THE UNIVERSITY OF CHICAGO 



SINGULARITIES, UNIVERSALITY, AND SCALING 
IN EVAPORATIVE DEPOSITION PATTERNS 



A DISSERTATION SUBMITTED TO 
THE FACULTY OF THE DIVISION OF THE PHYSICAL SCIENCES 
IN CANDIDACY FOR THE DEGREE OF 
DOCTOR OF PHILOSOPHY 

DEPARTMENT OF PHYSICS 



BY 

YURI O. POPOV 



CHICAGO, ILLINOIS 
DECEMBER 2003 



Copyright © 2008 by Yuri O. Popov 
All rights reserved 



To Denis 



ACKNOWLEDGMENTS 



I am greatly indebted to Tom Witten for his advice, everlasting support, and infinite 
patience. He has been an advisor, a teacher, and a friend during all these years; it is 
doubtful this work would be completed without him. I am deeply grateful. 

I am very thankful to Denis who was a source of inspiration during the latest years 
of this work, and who provided all the needed love and support. 

My parents gave me their trust and no-nonsense judgment that helped me to get 
to this stage; I would like to thank them. 

The University of Chicago was a wonderful place to work at during these years, 
and I would like to express my gratitude to many members of its community who 
made this experience truly unforgettable. Todd Dupont, David Grier, and Leo 
Kadanoff provided invaluable help during various stages of this work; other faculty 
members of the James Frank Institute contributed to the unique scientific atmo- 
sphere I was surrounded by. Postdocs and graduate students — Itai Cohen, Robert 
Deegan, Vladimir Belyi, Ajay Gopinathan, Toan Nguyen, Jung-Ren Huang, Ajay 
Gopal, Alexey Tkachenko, Joe Plewa, Brian DiDonna, Hiam Diamant, Swan Behrens, 
Matthias Mobius, Jennifer Curtis, Adrian Muresan, Dorel Buta, Li-Sheng Tseng, 
Shunji Egusa, Tao Liang — helped in a countless number of ways and generally made 
my experience at the university enjoyable, both academically and otherwise. 

My personal friends, both within and outside the university, — Bruno Carneiro 
da Cunha, Cristian Garcia, Pedro Amaral, Alexander Karaivanov, Mariana Gatzeva, 
Andrea Tiseno, Denis Suprun, Leonid Malyshkin, Andrey Tsvetkov, Andrey Tokarev, 
Denis Shpakov, Pavel Danilov, Boris Baryshnikov, Maria Razumova, Vladislav Gav- 
rilets — made this journey bright and exciting; my life without them would be only 
a fraction of the life I really had during these years. Thanks to all of them. 

Finally, I would like to thank Nobuko McNeill, Van Bistrow, Rosemary Garrison, 
and Melva Smith for easing both the formal and the informal parts of the university 
experience, for their care and indispensable assistance. 



iv 



TABLE OF CONTENTS 

ACKNOWLEDGMENTS iv 

LIST OF FIGURES vii 

ABSTRACT xi 

1 INTRODUCTION 1 

2 MAIN IDEAS AND GENERAL THEORY OF SOLUTE TRANSFER AND 
DEPOSIT GROWTH IN EVAPORATING DROPS 5 

2.1 System, geometry, coordinates, and assumptions 5 

2.2 On evaporation rate 7 

2.3 The full system of equations for hydrodynamic flow: conservation of 
mass, Darcy's law, and Young-Laplace equation 9 

2.4 Solute transfer and deposit growth 12 

3 DEPOSIT GROWTH FOR ZERO- VOLUME PARTICLES IN CIRCULAR 
EVAPORATING DROPS 15 

3.1 Surface shape 15 

3.2 Evaporation rate 16 

3.3 Flow field 22 

3.4 Solute transfer and deposit growth 23 

4 DEPOSIT GROWTH FOR ZERO- VOLUME PARTICLES IN POINTED EVAP- 
ORATING DROPS 25 

4.1 Surface shape 25 

4.2 Evaporation rate 28 

4.3 Flow field 32 

4.4 Streamlines 35 

4.5 Solute transfer and deposit growth 38 

5 DEPOSIT GROWTH FOR FINITE- VOLUME PARTICLES IN CIRCULAR 
EVAPORATING DROPS 49 

5.1 Model, assumptions, and geometry 49 

5.2 Main equations: conservation of mass 54 

5.3 Analytical results in the limit of small initial concentrations of the solute 58 

5.4 Numerical results for arbitrary initial concentrations of the solute . . 70 

5.5 Discussion 75 

6 CONCLUSIONS 79 



v 



vi 



A DOMINANCE OF THE CAPILLARY FORCES OVER VISCOUS STRESS, 
AND LEGITIMACY OF EMPLOYMENT OF THE EQUILIBRIUM SUR- 
FACE SHAPE 80 

REFERENCES 82 



LIST OF FIGURES 



2.1 Illustration of the possibility of a wide range of contact angles in the 
presence of a groove or another inhomogeneity 6 

2.2 The rate-limiting process for evaporative mass loss. It is the diffusion 
of saturated vapor just above the interface rather than the transfer 
across the interface 8 

2.3 Illustration of analogy between evaporation rate J for a liquid drop 
and electric field E for a conductor. Consideration of the drop (or 
conductor) and its reflection in the plane of the substrate significantly 
simplifies the boundary problem 9 

2.4 Conservation of mass: the liquid-vapor interface lowers exactly by the 
amount of fluid evaporated from the surface plus the difference between 

the outflow and the influx of fluid from the adjacent regions 10 

3.1 Numerical results: Dependence of contact angle 9 on time t. Different 
curves correspond to different initial contact angles; values of param- 
eter 9i are shown at each curve. The analytical result [Eq. (3.16)] in 
limit 9{ — > is also provided (the solid curve) 20 

3.2 Numerical results: Dependence of water mass M on time t. Different 
curves correspond to different initial contact angles; values of param- 
eter 9i are shown at each curve. The analytical result [Eq. (3.18)] in 
limit 9 1 — ■> is also provided (the solid curve) 20 

3.3 Dependence of the total drying time on the initial contact angle. The 
solid curve demonstrates the analytical result for small contact angles 
[Eq. (3.17)]; the dashed curve represents the exact numerical result 
[inferred from Eq. (3.15)] 21 

4.1 (a) A water droplet with a sector-shaped boundary on the plane sub- 
strate (side view), (b) The same droplet pictured from another point 
(top view). Black lines are the grooves on the substrate necessary to 
"pin" the contact line. (Courtesy Itai Cohen.) 26 

4.2 Geometry of the problem. The plane of the figure coincides with the 
substrate 27 

4.3 Dependence of exponent v in the power law h(r) [Eq. (4.2)] on opening 
angle a, after Ref. (26) 28 

4.4 Solid line: Dependence of exponent \x in the power law J(r) [Eqs. (4.5) 
and (4.12)] on opening angle a for an angular sector, after Refs. (36; 
37; 38). The same dependencies for a wedge (dashed line) and a cone 
(dotted line) are also shown 30 

4.5 Typical behavior of the numerical solution and the analytical asymp- 
totics of function ip((p) for two values of opening angle 34 



vii 



viii 



4.6 Characteristic behavior of flow field for two values of opening angle. 
Each arrow represents the absolute value and the direction of velocity 

v at the point of arrow origin 36 

4.7 Dependence of exponent e in the power law r(0) [Eq. (4.29)] on opening 
angle a. The two curves correspond to the two model forms for function 

J {(f) 38 

4.8 Typical shape of the streamlines for two values of opening angle (for 

the same values as in Figs. 4.5 and 4.6) 39 

4.9 Qualitative sketch: mutual location of streamlines (the two lines with 
arrows) and isochrones (the four numbered lines). Solute moves along 
the streamlines towards the contact line (the bold line). Shaded area 
is swept by an infinitesimal element of fluid between the two infinitesi- 
mally close streamlines as that element moves towards the contact line. 
The isochrones are the geometric locations, starting from which the so- 
lute reaches the contact line at the same time. Solute from isochrone 
1 reaches the contact line first; solute from isochrone 4 reaches the 
contact line last 41 

4.10 An illustration of the derivation of the four-thirds law at a straight 
segment of the contact line. The contact line is normal to the plane of 
the figure. Length drQ is along the contact line and hence not shown. 

The flow is in the plane of the figure from left to right 43 

4.11 Exponent of distance tq in the power law dm/dro(ro) [Eqs. (4.37) and 
(4.40)] as a function of the opening angle for the two time regimes. 
The early-time curve corresponds to the exponent (3 of Eq. (4.38); the 
intermediate-time curves correspond to the exponent 7 of Eq. (4.42). 
The two curves for the intermediate-time exponent correspond to the 

two model forms for function J(4>) 45 

4.12 Exponent of time t in the power law dm/drQ(t) [Eqs. (4.37) and (4.40)] 
as a function of the opening angle for the two time regimes. The 
early-time curve corresponds to the exponent 2/(1 + A) = 4/3; the 
intermediate-time curves correspond to the exponent 5 of Eq. (4.41). 
The two curves for the intermediate-time exponent correspond to the 

two model forms for function J(0) 46 

4.13 Log-log plot of the numerical solution d 2 m/ 'dtdr§{t) for two values of 
opening angle (for the same values as in Figs. 4.5, 4.6 and 4.8) 47 

4.14 Log-log plot of the numerical solution d^m/ 'dtdr-Q^o) for two values of 
opening angle (for the same values as in Figs. 4.5, 4.6, 4.8 and 4.13). 
Parameter r* is defined by Eq. (4.44) 48 

5.1 Mutual location of two "phases" in the drying drop: L is the liquid 

phase, and D is the deposit phase 51 



ix 



5.2 Time evolution of the deposit phase growth: side view (left) and top 
view (right). Only the deposit phase is shown. Thickness of the ring 

is exaggerated compared to the typical experimental results 51 

5.3 Presence of particles in the deposit does not obstruct fluid evaporation 
at the edge of the drop. All the necessary fluid is supplied, and it is 
this motion of the fluid that brings the particles to the deposit phase. 
Also shown schematically is the fact that the boundary between the 
phases is vertical and the particles get stacked at full height between 

the substrate and the free surface of the drop 52 

5.4 Geometry of the problem. Vertical scale is exaggerated in order to see 

the details, typically H <C and h <C 53 

5.5 Mass of a drying drop as a function of time. Experimental results, 
after Refs. (21; 22). The line running through the data is a linear fit. 
(Courtesy Robert Deegan.) 62 

5.6 Fraction of the solute in the deposit phase v£/V S as a function of 
time [Eq. (5.49)] 63 

5.7 Height of the phase boundary H(t) in units of 8^R^ \Jxilv as a function 

of time [Eq. (5.50)] 64 



5.8 Width of the deposit ring W(t) in units of RiyXilP as a function of 
time [Eq. (5.51)] 65 

5.9 Deposit ring profile: height vs. width. The vertical scale is different 
from the horizontal scale by a factor of 8^ <C 1 65 

5.10 Ring width vs. drop radius. Experimental results, after Refs. (21; 22). 
The line running through the data is a linear fit. (Courtesy Robert 
Deegan.) 67 

5.11 Ring width normalized by the drop radius vs. initial concentration of 
the solute for two different particle sizes. Experimental results, after 
Refs. (21; 22). The two data sets are offset by a factor of 5 to avoid 
mixing of the data points related to the different particle sizes. The 
lines running through the data are linear fits in the double-logarithmic 
scale, which upon conversion to the linear scale yield power laws with 
exponents 0.78 ± 0.10 and 0.86 ± 0.10. (Courtesy Robert Deegan.) . . 68 

5.12 A photographic sequence demonstrating a depinning event. Experi- 
mental results, after Refs. (21; 22). The view is from above, and the 
solid white band in the lower part of the frame is the ring; the rest 
of the drop is above the ring. The time between the first and the last 
frames is approximately 6 s; the major axis of the hole is approximately 

150 jum. (Courtesy Robert Deegan.) 68 



5.13 Depinning time normalized by the extrapolated drying time vs. ini- 
tial concentration of the solute. Experimental results, after Refs. (21; 
22). The line running through the data is a linear fit in the double- 
logarithmic scale, which upon conversion to the linear scale yields a 
power law with exponent 0.26 ± 0.08. (Courtesy Robert Deegan.) . . 69 

5.14 Numerical results: Dependence of angle between the liquid-air in- 
terface and the substrate at the phase boundary on time t. Different 
curves correspond to different initial concentrations of the solute; val- 
ues of parameter Xi/p are shown at each curve. The analytical result 

in limit Xi/p ~^ is a lso provided 71 

5.15 Numerical results: Dependence of the volume fraction of the solute in 
the deposit phase Vf^/V^ on time t. Different curves correspond to 
different initial concentrations of the solute; values of parameter Xi/p 
are shown at each curve. The analytical result in limit Xi/p ~^ is 

also provided 72 

5.16 Numerical results: Dependence of the height of the phase boundary H 
on time t. Different curves correspond to different initial concentrations 
of the solute; values of parameter Xi/p are shown at each curve. The 
analytical result in limit Xi/p ~^ is also provided 73 

5.17 Numerical results: Dependence of the width of the deposit ring W on 
time t. Different curves correspond to different initial concentrations 
of the solute; values of parameter Xi/p are shown at each curve. The 
analytical result in limit Xi/p ~ ¥ is a lso provided 73 

5.18 Numerical results: Ring profile. Dependence of the height of the phase 
boundary H on the width of the deposit ring W . Different curves 
correspond to different initial concentrations of the solute; values of 
parameter Xi/p are shown at each curve. The analytical result in limit 
Xi/p is also provided 74 

5.19 Numerical results: Log-log plot of the dependence of the height of 
the phase boundary H and the width of the deposit ring W on the 
initial volume fraction of the solute. "Lg" stands for the logarithm of 
base 10, and thus the numbers along the horizontal axis represent the 
order of magnitude of the ratio Xi/P- The main-order analytical results 

H oc J Xi/p an d W oc J Xi/p are a lso provided for comparison 74 



ABSTRACT 



The theory of solute transfer and deposit growth in evaporating sessile drops on a 
plane substrate is presented. The main ideas and the principal equations are formu- 
lated. The problem is solved analytically for two important geometries: round drops 
(drops with circular boundary) and pointed drops (drops with angular boundary). 
The surface shape, the evaporation rate, the flow field, and the mass of the solute 
deposited at the drop perimeter are obtained as functions of the drying time and the 
drop geometry. In addition, a model accounting for the spatial extent of the deposit 
arising from the non-zero volume of the solute particles is solved for round drops. 
The geometrical characteristics of the deposition patterns as functions of the drying 
time, the drop geometry, and the initial concentration of the solute are found ana- 
lytically for small initial concentrations of solute and numerically for arbitrary initial 
concentrations of solute. The universality of the theoretical results is emphasized, 
and comparison to the experimental data is made. 
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CHAPTER 1 
INTRODUCTION 



The problem of the so-called "coffee-drop deposit" has recently aroused great inter- 
est. The residue left when coffee dries on the countertop is usually darkest and hence 
most concentrated along the perimeter of the stain. Ring-like stains, with the so- 
lute segregated to the edge of a drying drop, are not particular to coffee. Mineral 
rings left on washed glassware, banded deposits of salt on the sidewalk during win- 
ter, and enhanced edges in water color paintings are all examples of the variety of 
physical systems displaying similar behavior and understood by coffee-drop deposit 
terminology. 

Understanding the process of drying of such solutions is important for many sci- 
entific and industrial applications, where ability to control the distribution of the 
solute during drying process is at stake. For instance, in the paint industry, the pig- 
ment should be evenly dispersed after drying, and the segregation effects are highly 
undesirable. Also, in the protein crystallography (1; 2), attempts are made to assem- 
ble the two-dimensional crystals by using evaporation driven convection, and hence 
solute concentration gradients should be avoided. On the other hand, in the pro- 
duction of nanowires (3) or in surface patterning (4) perimeter-concentrated deposits 
may be of advantage. Recent important applications of this phenomenon related to 
DNA stretching in a flow have emerged as well (5). For instance, a high-throughput 
automatic DNA mapping was suggested (6), where fluid flow induced by evaporation 
is used for both stretching DNA molecules and depositing them onto a substrate. 
Droplet drying is also important in the attempts to create arrays of DNA spots for 
gene expression analysis. 

Ring-like deposit patterns have been studied experimentally by a number of groups. 
Difficulties of obtaining a uniform deposit (7), deformation of sessile drops due to a 
sol-gel transition of the solute at the contact line (8; 9), stick-slip motion of the 
contact line of colloidal liquids (10; 11), and the effect of ring formation on the evap- 
oration of the sessile drops (12) were all reported. The evaporation of the sessile drops 
(regardless of solute presence) has also been investigated extensively. Constancy of 
the evaporation flux was demonstrated (13; 14), and the change of the geometrical 
characteristics (contact angle, drop height, contact-line radius) during drying was 
measured in detail (15; 16; 17; 18). 

The most recent and complete experimental effort to date on coffee-drop deposits 
was conducted by Robert Deegan et al. (19; 20; 21; 22). Most experimental data 
referred to in this work originate from observations and measurements of this group. 
They reported extensive results on ring formation and demonstrated that these could 
be quantitatively accounted for. The main ideas of the theory of solute transfer in 
such physical systems have also been developed in their work (19). It was observed 
that the contact line of a drop of liquid remains pinned during most of the drying 
process. While the highest evaporation occurs at the edges, the bulk of the solvent is 
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concentrated closer to the center of the drop. In order to replenish the liquid removed 
by evaporation at the edge, a flow from the inner to the outer regions must exist inside 
the drop. This flow is capable of transferring all of the solute to the contact line 
and thus accounts for the strong contact-line concentration of the residue left after 
complete drying. The idea of this theory is very robust since it is independent of the 
nature of the solute and only requires the pinning of the edge during drying (which can 
occur in a number of possible ways: surface roughness, chemical heterogeneities etc). 
This theory accounts quantitatively for many phenomena observed experimentally; 
among other things, we will reproduce its main conclusions in this work. 

Mathematically, the most complicated part of this problem is related to determin- 
ing the evaporation rate from the surface of the drop. An analogy between the diffu- 
sive concentration fields and the electrostatic potential fields was suggested (23; 24), 
so that an equivalent electrostatic problem can be solved instead of the problem of 
evaporation of a sessile drop. Important analytical solutions to this equivalent elec- 
trostatic problem in various geometries were first derived by Lebedev (23). A number 
of useful consequences from these analytical results were later reported in Ref. (25). 

In this work, we discuss the theory of solute transfer and deposit growth in evap- 
orating sessile droplets on a substrate and provide quantitative account for many 
observed phenomena and measurement results. Chapter 2 discusses the main ideas 
and the general theory; all principal equations are derived in that chapter. While 
most of its equations have been reported previously, the derivation presented here is 
original and deals with some mathematical issues never fully addressed before in the 
context of the current problem. Chapter 2 is the basis for all the following chapters 
of this work, and its content is required for all the other chapters. 

While the principal equations are fully derived and presented in Chapter 2, their 
solution depends heavily on the geometry of the drop. The flow pattern discussed in 
this work is a type of hydrodynamic flow that is sensitive to the perimeter shape, i.e. 
the shape of the contact line. Mathematically, solution to the differential equations 
depends on the boundary conditions. Chapters 3 and 4 discuss the analytical solution 
to this problem in two important geometries. The two geometries are the drops 
with circular boundary (round drops) and the drops with angular boundary (pointed 
drops). The choice of these two geometries is not accidental. An arbitrary boundary 
line can be represented as a sequence of smooth segments, which can be approximated 
by circular arcs, and fractures, which can be approximated by angular regions. Thus, 
knowledge of analytical solution for both circular (Chapter 3) and angular (Chapter 4) 
boundary shapes fills out the quantitative picture of solute transfer and deposit growth 
for an arbitrary drop boundary. 

The case of the round drops is the most important from the practical point of 
view and the easiest to deal with mathematically. This case allows for a full analytical 
solution, and this solution has been obtained earlier. Here, in Chapter 3, we reproduce 
concisely the earlier results and report some new ones. Its content is a prerequisite 
to Chapter 5, but is not essential for understanding Chapter 4 (although it is used 
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for drawing some parallels and for comparison of the results in the two geometries). 

The case of the pointed drops (Chapter 4), while also important, is much more 
complicated mathematically than the round-drop solution is. Presence of the vertex 
of the angle introduces a singularity at this vertex in addition to the weaker singularity 
at the contact line. Singularities govern the solutions to differential equations, and 
thus presence of the angle and its vertex changes the results substantially. Also, an 
angular region, as a mathematical object, is infinite, while a circular region is always 
bounded. The real drops with a fracture must always have a third (the furthest from 
the vertex) side of its contact line, and therefore the overall solution depends on the 
shape of that furthest part of the boundary. At the same time, we are interested 
only in the universal features of the solution that are independent of that furthest 
part. Keeping in mind this lust for universality and the mathematical complexity 
mentioned above, we specify only one boundary of the drop (the vertex and the two 
sides of the angle) leaving the remainder of the boundary curve unspecified. Such an 
approach turns out to be sufficient to determine the universal features of the solution, 
and it allows us to find all the important singularities as power laws of distance from 
the vertex of the angle. Most of the results of Chapter 4 were originally obtained 
in our earlier works (26; 27). Chapters 4 and 5 are completely independent of each 
other. 

Chapters 2-4 address the issue of the deposit mass accumulation at the drop 
boundary, however, they treat the solute particles as if they do not occupy any volume, 
and hence all the solute can be accommodated at the one-dimensional singularity of 
the contact line. In reality, the solute deposit accumulated at the perimeter of a 
drying drop has some thickness, for instance, the shape of the solute residue in a 
round drop is more like a ring rather than an infinitely thin circumference of the 
circle. We attribute this finite volume of the solute deposit simply to the finite size of 
the solute particles, i.e. we assume the particles do occupy some volume and hence 
cannot be packed denser than certain mass per unit volume. A model accounting 
for the finite size of the deposit and determining its geometric characteristics (height, 
width) is considered in Chapter 5. The model is solved for the simplest case of circular 
geometry, and its results are compared to the zero-volume results of Chapter 3 and 
to the experimental results of Refs. (21; 22). Both the analytical and the numerical 
solutions are provided, and both compare well with the experimental data (and with 
each other). Results of Chapter 5 are presented for the first time. A chapter of 
conclusions completes this work. 

Before proceeding to the main matter, we would like to point out some features 
common to all our results and giving rise to the title of this work. Flows found here 
are capable of transferring 100% of the solute to the contact line, and thus account for 
the strong perimeter concentration of the solute in all cases. Many quantities scale as 
power laws; some others follow different functional dependencies. However, both the 
exponents of the power laws and any elements of the other functional dependencies 
turn out to be universal within our model. They do not depend on any parameters of 
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the system other than the system geometry. Within the range of applicability of our 
theory, there are no fitting parameters, no undetermined constants, and no unknown 
coefficients. Thus, for instance, exponents of the power laws are as universal as the 
exponent of distance —2 in the Coulomb's law. We view this universality as one of 
the main advantages and one of the most exciting features of this theory. 



CHAPTER 2 

MAIN IDEAS AND GENERAL THEORY OF SOLUTE 
TRANSFER AND DEPOSIT GROWTH IN 
EVAPORATING DROPS 

In this chapter we will present the general theory of solute transfer to the contact line 
that will subsequently be solved analytically for two important geometries. 

2.1 System, geometry, coordinates, and assumptions 

We consider a sessile droplet of solution on a horizontal surface (substrate). The 
nature of the solute is not essential for the mechanism. The typical diameter of the 
solute particles in Deegan's experiments was of the order of 0.1-1 /im; we will assume 
a similar order of magnitude throughout this work. For smaller particles diffusion 
becomes more important compared to the hydrodynamic flows of this work. For 
larger particles sedimentation may become important when particles exceed certain 
size. 

The droplet is bounded by the contact line in the plane of the substrate. The 
(macroscopic) contact line is defined as the common one-dimensional boundary of 
all three phases (liquid, air and solid substrate). We will not specify the shape of 
the contact line in this chapter; in the later chapters we will make this selection and 
describe the two qualitatively different limits — circular drops and angular drops. 
As we explained in the introduction, these shapes account for the smooth and frac- 
tured segments of an arbitrary contact line. Main equations are not sensitive to the 
particular geometry, only the boundary conditions are. 

We assume that the droplet is sufficiently small so that the surface tension is 
dominant, and the gravitational effects can be safely neglected. Mathematically, this 
is controlled by the Bond number Bo = pgRh max /cr, which accounts for the balance 
of surface tension and gravitational force on the surface shape. Here p is the density 
of fluid, g is the gravitational constant, R is a typical size of the drop in the plane 
of the substrate, h max is the maximal height of the drop at the beginning of the 
drying process, and a is the surface tension at the liquid-air interface. For the typical 
experimental conditions the Bond number is of the order of 0.02-0.05, and thus gravity 
indeed is unimportant and surface shape is governed mostly by the surface tension. 

At the same time, we do not assume that the contact angle between the liquid- 
air interface and the plane of the substrate is constant along the contact line on 
the substrate, nor do we assume it is constant in time. To achieve a prescribed 
boundary shape (other than a perfect circle on an ideal plane), the substrate must 
have scratches, grooves or other inhomogeneities (sufficiently small compared to the 
dimensions of the droplet), which pin the contact line. A strongly pinned contact line 
can sustain a wide range of (macroscopic) contact angles. The contact angle is not 
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Substrate 



Figure 2.1: Illustration of the possibility of a wide range of contact angles in the 
presence of a groove or another inhomogeneity. 



fixed by the interfacial tensions as it is on a uniform surface (Fig. 2.1). Throughout 
this work we will deal with small contact angles (9 <^ 1) as is almost always the case 
in the experimental realizations (typically, max < 0.1-0.3); however, the general 
equations of this chapter do not rely on the smallness of the contact angle. 

We will use the cylindrical coordinates (r, 0, z) throughout this work, as most 
natural for both geometries of interest. Coordinate z is always normal to the plane of 
the substrate, and the plane itself is described by z = 0, with z being positive on the 
droplet side of the space. Coordinates (r, 0) are the polar radius and the azimuthal 
angle, respectively, in the plane of the substrate. The origin is chosen at the center 
of the circle in the circular geometry and at the vertex of the angle in the angular 
geometry. The geometry of the problem is quite complicated despite the visible 
simplicity. In both geometries of interest, we consider an object, whose symmetry 
does not match the symmetry of any simple orthogonal coordinate system of the three- 
dimensional space. For instance, solution of the Laplace equation (needed below) 
requires introduction of the special coordinate systems (the toroidal coordinates and 
the conical coordinates) with heavy use of various special functions. Similar difficulties 
related to the geometry arise in the other parts of the problem as well. 

We describe the surface shape of the drop h(r, <fi) by local mean curvature K 
that is spatially uniform at any given moment of time, but changes with time as 
droplet dries. Ideally, the surface shape should be considered dynamically together 
with the flow field inside the drop. However, as we show below, for flow velocities 
much lower than the characteristic velocity v* = a/3r] (where a is surface tension and 
i] is dynamic viscosity), which is about 24 m/s for water under normal conditions, 
one can consider the surface shape independently of the flow and use the equilibrium 
result at any given moment of time for finding the flow at that time. Another way of 
expressing the same condition is to refer to the capillary number Ca = rjv/a (where 
v is the characteristic value of the flow velocity of the order of 1-10 jLtm/s), which is 
the ratio of viscous to capillary forces. This ratio is of the order of 10~ 8 -10~ 7 under 
typical experimental conditions, clearly demonstrating that capillary forces are by far 
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the dominant ones in this system and that surface shape is practically equilibrium 
and depends on time adiabatically. 

We consider slow flows, i.e. flows with low Reynolds numbers (also known as 
"creeping flows"). This amounts to the neglect of the inertial terms in the Navier- 
Stokes equation. As all the conditions above are, this is well justified by the real 
experimental conditions. We also employ the so-called "lubrication approximation" . 
It is essentially based on the two conditions reflecting the thinness of the drop and 
resulting from the separation of the vertical and horizontal scales. One is that the 
pressure inside the drop p does not depend on the coordinate z normal to the sub- 
strate: d z p = 0. The other is related to the small slope of the free surface \Vh\ <C 1, 
which is equivalent to the dominance of the ^-derivatives of any component of flow 
velocity u: d z Ui ^> d s Ui (index s refers to the derivatives with respect to any coor- 
dinate in the plane of the substrate). The lubrication approximation is a standard 
simplifying procedure for this class of hydrodynamic problems (28; 29; 30). 

Before proceeding to the main section of this chapter and formulating the main 
ideas of the theory, we will make a brief note on evaporation rate. 

2.2 On evaporation rate 

In order to determine the flow caused by evaporation, one needs to know the flux 
profile of liquid leaving each point of the surface by evaporation. This quantity will 
be seen to be independent of the processes going on inside the drop and must be 
determined prior to considering any such processes. 

The functional form of the evaporation rate J(r, 0) (defined as evaporative mass 
loss per unit surface area per unit time) depends on the rate-limiting step, which 
can, in principle, be either the transfer rate across the liquid-vapor interface or the 
diffusive relaxation of the saturated vapor layer immediately above the drop. We 
assume everywhere that the rate-limiting step is diffusion of liquid vapor (Fig. 2.2) 
and that evaporation rapidly attains a steady state. Indeed, the transfer rate across 
the liquid-vapor interface is characterized by the time scale of the order of 10~ 10 s, 
while the diffusion process has characteristic times of the order of R^/D (where D 
is the diffusion constant for vapor in air and R is a characteristic size of the drop), 
which is of the order of seconds for water drops under typical drying conditions. Also, 
the ratio of the time required for the vapor-phase water concentration to adjust to 
the changes in the droplet shape (R 2 /D) to the droplet evaporation time tf is of the 
order of (n s — noo)/ p ~ 1CP 5 , where n s is the density of saturated vapor just above 
the liquid-air interface and is the ambient vapor density (25). Thus, indeed, vapor 
concentration adjusts rapidly compared to the time required for water evaporation, 
and the evaporation process can be considered quasi-steady. 

As the rate-limiting process is the diffusion, the density of vapor above the liquid- 
vapor interface n obeys the diffusion equation. Since diffusion rapidly attains a steady 
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Figure 2.2: The rate-limiting process for evaporative mass loss. It is the diffusion of 
saturated vapor just above the interface rather than the transfer across the interface. 



state, this diffusion equation reduces to the Laplace equation 

V 2 n = 0. (2.1) 

This equation is to be solved together with the following boundary conditions de- 
pendent on the geometry of the drop: (a) along the surface of the drop the air is 
saturated with vapor and hence n at the interface is the constant density of saturated 
vapor n s , (b) far away from the drop the density approaches the constant ambient 
vapor density rico, and (c) vapor cannot penetrate the substrate and hence d z n = 
at the substrate outside of the drop. Having found density of vapor, one can obtain 
the evaporation rate J = —DWn, where D is the diffusion constant. 

This boundary problem is mathematically equivalent to that of a charged conduc- 
tor of the same geometry at constant potential if we identify n with the electrostatic 
potential and J with the electric field. Moreover, since there is no component of J 
normal to the substrate, we can further simplify the boundary problem by considering 
a conductor of the shape of our drop plus its reflection in the plane of the substrate 
in the full space instead of viewing only the semi-infinite space bounded by the sub- 
strate (Fig. 2.3). This reduces the number of boundary conditions to only two: (a) 
n = n s on the surface of the conductor, and (b) n = rioo at infinity. The shape of the 
conductor (the drop and its reflection in the substrate) is now symmetric with respect 
to the plane of the substrate. This plane-symmetric problem of finding the electric 
field around the conductor at constant potential in infinite space is much simpler than 
the original problem in the semi-infinite space, and this will be the problem we will 
actually be solving in order to find the evaporation rate for the two geometries of 
interest. Particular solution in each of the geometries will be presented in the next 
two chapters. 

Having formulated physical assumptions intrinsic to the theory, we are now in 
position to state its main ideas. 
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Figure 2.3: Illustration of analogy between evaporation rate J for a liquid drop and 
electric field E for a conductor. Consideration of the drop (or conductor) and its 
reflection in the plane of the substrate significantly simplifies the boundary problem. 

2.3 The full system of equations for hydrodynamic flow: 
conservation of mass, Darcy's law, and Young-Laplace 

equation 

The essential idea behind the theory has been developed in works of Deegan et al. (19; 
20; 21). It is an experimental observation that the contact line of a drop of liquid is 
pinned during most of the drying process. While the highest evaporation occurs at 
the edges, the bulk of the solvent is concentrated closer to the center of the drop. In 
order to replenish the liquid removed by evaporation at the edge, a flow from the inner 
to the outer regions must exist inside the drop. This flow is capable of transferring 
all of the solute to the contact line and thus accounts for the strong contact-line 
concentration of the residue left after complete drying. Thus, a pinned contact line 
entails fluid flow toward that contact line. The "elasticity" of the liquid-air interface 
fixed at the contact line provides the force driving this flow. 

To develop this idea mathematically, we ignore for a moment any solute in the 
liquid. Once the flow is found, one can track the motion of the suspended particles, 
since they are just carried along by the flow. The purpose of this section is to describe 
a generic method for finding the hydrodynamic flow of the liquid inside the drop, 
which is a prerequisite to knowing the details of solute transfer. 

We define depth-averaged flow velocity by 



where u s is the in-plane component of the local three-dimensional velocity u. Then 
we write the conservation of fluid mass in the form 




(2.2) 




(2.3) 




Figure 2.4: Conservation of mass: the liquid-vapor interface lowers exactly by the 
amount of fluid evaporated from the surface plus the difference between the outflow 
and the influx of fluid from the adjacent regions. 



where t is the time, p is the density of the fluid, and each of the quantities h, J and v is 
a function of r, and t. (We will drop the (V/i) 2 part of the second term everywhere 
in the following since it is always small compared to unity, as will be seen below.) 
This equation represents the fact that the rate of change of the amount of fluid in a 
volume element (column) above an infinitesimal area on the substrate (third term) is 
equal to the negative of the sum of the net flux of liquid out of the column (first term) 
and the amount of mass evaporated from the surface element on top of that column 
(second term); Fig. 2.4 illustrates this idea. Thus, this expression relates the depth- 
averaged velocity field v(r, 0, t) to the liquid-vapor interface position h(r, 0, t) and 
the evaporation rate J(r, 0, t) . However, this is only one equation for two variables 
since vector v has generally two components in the plane of the substrate. Moreover, 
while the evaporation rate J is indeed independent of flow v, the free-surface shape 
h should in general be determined simultaneously with v. Thus, there are actually 
three unknowns to be determined together (h and two components of velocity, say, v r 
and va), and hence two more equations are needed. These additional equations will 
be of the hydrodynamic origin. 

We start with the Navier-Stokes equation with inertial terms omitted (low Rey- 
nolds numbers): 

Vp = ??V 2 u, (2.4) 

where p is the fluid pressure, 77 is the dynamic viscosity, and u is the velocity. Applying 
lubrication- approximation conditions d z p = and d z U{ ^> d a u^ we arrive at the 
simplified form of this equation 

V s p = r]d zz u s , (2.5) 
where index s again refers to the vector components along the substrate. From now on 
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we will suppress the subscript s at the symbol of nabla-operator, and will assume for 
the rest of this work that this operator refers to the two-dimensional vector operations 
in the plane of the substrate. Solution to the above equation with boundary conditions 

u s | z=0 = and d z u a \ z=h = (2.6) 

(no slip at the substrate and no stress at the liquid-air interface) yields 



or, after vertical averaging (2.2), 

h 2 

v = -f Vp. (2.8) 

This result is a variant of the Darcy's law (30; 31). Note that since a curl of a gradient 
is always zero and rj is a constant, the preceding equation can be re- written as 

Vx(J)= 0. (2.9) 

This condition is analogous to the condition of the potential flow (V x v = 0), but 
with a quite unusual combination of the velocity and the surface height v//i 2 in place 
of the usual velocity v. 

Relation (2.8) provides the two sought equations in addition to the conservation 
of mass (2.3). However, it contains one new variable, pressure p, and hence another 
equation is needed. This last equation is provided by the condition of the mechanical 
equilibrium of the liquid-air interface (also known as the Young-Laplace equation) 
relating the pressure and the surface shape: 

p=-2Ka + p atm . (2.10) 

Here p a tm is the atmospheric pressure, a is the surface tension, and K is the mean 
curvature of the surface, uniquely related to the surface shape h by differential geom- 
etry. Note that this expression is independent of both the conservation of mass (2.3) 
and the Darcy's law (2.8). Thus, the complete set of equations required to fully 
determine the four dynamic variables h, p, v r , and iu consists of four differential 
equations (together with the appropriate boundary conditions at the contact line, 
which are dependent on the particular geometry of the drop): one equation of the 
conservation of mass (2.3), two equations of the Darcy's law (2.8), and one equation 
of the mechanical equilibrium of the interface (2.10). They provide all the necessary 
conditions to solve the problem at least in principle. 

In practice, however, solution of these four coupled differential equations is not 
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possible in most geometries of practical interest. At the same time, under normal 
drying conditions the viscous stress is negligible, or, equivalently, the typical velocities 
are much smaller than v* = u/Stj ~ 24 m/s (for water under normal conditions). As 
we show in the Appendix, the four equations decouple under these conditions. As 
a result, one can employ the equilibrium result for the surface shape h at any given 
moment of time, and then determine the pressure and the velocity fields for this 
fixed functional form of h. Mathematically, the original system of equations can be 
rewritten as: 

2K = -^, (2.11) 
o~ 

V • (h 3 Vij) = --- d t h, (2.12) 

p 

v = h 2 Vij, (2.13) 

where Ap = p^—Patmi = ~ e Pl/3T], and p$ and p\ are the leading and the first-order 
terms in the expansion of pressure p = Po + ep\ + • • • in a small parameter e inversely 
proportional to v* (see the Appendix for details). Note that pq is independent of (r, 0), 
although it does depend on time (this time dependence will be determined later in this 
work). Therefore, there is a profound difference between equations (2.10) and (2.11): 
the former is a local statement, with the right-hand side depending on the coordinates 
of a point within the drop, while the latter is a global condition of spatial constancy 
of the mean curvature throughout the drop. Equation (2.11) defines the equilibrium 
surface shape for any given value of pq at any given moment of time, and moreover, 
can be solved independently of the other equations. Thus, the procedure for finding 
the solution becomes significantly simplified: first find the equilibrium surface shape h 
from condition (2.11) and independently specify the functional form of the evaporation 
rate J from an equivalent electrostatic problem, then solve equation (2.12) for the 
reduced pressure ip, and finally obtain the flow field v according to prescription (2.13). 
The next two chapters will be devoted to the particular steps of this procedure for 
the two geometries of interest. 

In the next section, we will describe how knowledge of flow inside the drop allows 
one to find the rate of solute transfer to the contact line and determine the laws of 
deposit growth. 

2.4 Solute transfer and deposit growth 

With the velocity field inside the drop in hand, we can compute the rate of the deposit 
growth at the contact line. We assume that the suspended particles are carried along 
by the flow with velocity equal to the fluid velocity. Integrating the velocity field: 

dr _ v r (r,(f)) ^ 



rd(p ^(r,0)' 
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we find the streamline equation r(fa) or far), i.e. the trajectory of each particle 
as it moves with the fluid. This streamline equation is independent of the overall 
intensity of evaporation [since both v r and depend on it as a multiplicative factor 
dropping out of Eq. (2.14)], and thus the shape of the streamlines is universal for each 
geometry of the drop. Physically, this indicates that solute particles move along the 
same trajectories independently of how fast evaporation occurs and hence how fast 
the flow is. 

Given the shape of the streamlines, we can compute the time it takes an element 
of fluid (having started from some initial point (r^, fa) and moving along a streamline) 
to reach the contact line. This time can be found by integrating both sides of either 
v r (r,<f),t)dt = dr or vAr,(f),t)dt = rd<p with known dependences of v r or va on the 
variables and known relation between r and <f) on the streamline. The integrations 
are to be conducted from r\ or fa (at time 0) to rj or 0^ (at time t), respectively, 
where (rf,(f)f) is the terminal endpoint of the trajectory on the contact line. Thus, 
the initial location of particles that reach the contact line at time t is characterized by 
(rj, fa). First, only particles initially located near the contact line reach that contact 
line. As time goes by, particles initially located further away from the contact line 
and in the inner parts of the drop reach the contact line. Finally, particles initially 
located in the innermost parts of the drop (e.g. at the center of a round drop or near 
the bisector of an angular drop) reach the contact line as well. The more time elapsed, 
the more particles reached the contact line and the larger the area is where they were 
spread around initially. One can view this process as inward propagation of the inner 
boundary of the set of initial locations of the particles that have reached the contact 
line by time t. As is easy to understand, the velocity of this front is equal to the 
negative of the vector of fluid velocity at each point (the fluid and the particles move 
towards the contact line while this front moves away from it, hence a minus sign). 

Within the time computed in the preceding paragraph all the solute that lays 
on the way of an element of fluid as it moves toward the contact line becomes part 
of the deposit. Now, we use our knowledge of the initial distribution of the solute, 
namely, that the solute has constant concentration c everywhere in the drop at time 
t — 0, and compute the mass of the deposit accumulated at the contact line by that 
time. The mass of the deposit dm accumulated at the contact-line element of length 
dl can be found by integrating h(r, 0, 0) over area dA between two infinitesimally 
close streamlines (terminating dl apart on the contact line) swept by the infinitesimal 
element of fluid on its way to the contact line and multiplying the result by the initial 
concentration c of the solute: 

dm = cf h(r,faO)rdrd<f>. (2.15) 
J dA 

Obviously, this mass will depend of the initial location (r^, fa) of the element of fluid. 
Thus, both the time elapsed from the beginning of the drying process (t) and the 
deposit mass accumulated at the contact line (dm) depend on the initial coordinates 
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of the arriving element of fluid (only one of the coordinates is actually an independent 
variable — the other is constrained by the streamline equation). Eliminating these 
coordinates from the expressions for time and mass, one can finally obtain the deposit 
mass accumulated at the contact line as a function of time elapsed from the beginning 
of the drying process. 

Thus, the procedure described allows one to find the rate of solute transfer and 
the laws of deposit growth. Since we use depth-averaged velocity throughout this 
work, we implicitly assume that there is no vertical segregation of the solute. 



CHAPTER 3 

DEPOSIT GROWTH FOR ZERO-VOLUME PARTICLES 
IN CIRCULAR EVAPORATING DROPS 

In this chapter we will provide the full solution to the problem of solute transfer 
in the case of round drops. In this geometry, the contact line is the circumference 
of a circle (and the origin of the cylindrical coordinates is located at the center of 
that circle). This case is of most importance from the point of view of practical 
applications and at the same time is the easiest to treat analytically. Most results 
can be derived in a closed analytical form, and many of them have been obtained in 
earlier works (19; 25). Some of them, however, are presented for the first time, and 
some correct earlier expressions. 



3.1 Surface shape 

We will follow the procedure explained in great detail in the preceding chapter. First, 
one needs to determine the equilibrium shape of the liquid-air interface. In circular 
geometry this task is particularly easy. Let R be the radius of the circular projection 
of the drop onto the plane of the substrate and 9 be the (macroscopic) contact angle. 
Then the solution to Eq. (2.11) with boundary condition h(R,(f),t) = is just a 
spherical cap. In cylindrical coordinates function h(r, 0, t) is independent of and 
can be written as 

h(r,t) = -r 2 -Root 6, (3.1) 

V shr 9 

where the radius of the footprint of this cap R and the contact angle 9 are related 
via the right-hand side of Eq. (2.11): 

R=^L s i n 0. (3.2) 
Ap 

Both 6 and a/Ap change with time during drying process; however, the radius of 
the footprint R stays constant. In most experimental realizations < 1 and the 
preceding expression takes even simpler form: 

2 2 

h(r,t) = ^-9(t)+0(9 3 ). (3.3) 

Thus, for small 9, quantity Vh is indeed small (since r < R) and can be safely 
neglected with respect to unity {i.e. the free surface of the drop is nearly horizontal) 
as was asserted in the preceding chapter. 

Knowledge of the surface shape allows one to find all the necessary geometrical 
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characteristics of the drop, for instance, its volume 

V = J* h{r, t) 1-nrdr = ttR 3 V{9) (3.4) 

or the total mass of the water (or any other fluid the drop is comprised of) 

M = pV = irpR 3 V(6), (3.5) 

where p is the density of the water and V(9) is a function of the contact angle: 

cos3* 3cos* + 2 * 3) m 
3 sm° V 4 

(the last equality is an expansion in limit 6 <C 1). 

3.2 Evaporation rate 

The other prerequisite to determining the flow inside the drop is the evaporation 
rate from the free surface of the drop. This task involves solution of the equivalent 
electrostatic problem (the Laplace equation) for the conductor of the shape of the 
drop plus its reflection in the plane of the substrate (kept at constant potential, as 
a boundary condition). In the case of the round drop the shape of this conductor 
resembles a symmetrical double-convex lens comprised of two spherical caps. The 
system of orthogonal coordinates that matches the symmetry of this object (so that 
one of the coordinate surfaces coincides with the surface of the lens) is called the 
toroidal coordinates (a, j3, 0), where coordinates a and (5 are related to the cylindrical 
coordinates r and z by 

R sinh a Rsinf3 ^ ^ 



cosh a — cos (5 cosh a — cos (5 

and the azimuthal angle has the same meaning as in cylindrical coordinates. Solu- 
tion to the Laplace equation in toroidal coordinates involves the Legendre functions 
of fractional degree and was derived in a book by Lebedev (23). The expression for 
the electrostatic potential or vapor density in toroidal coordinates obtained in that 
book is independent of the azimuthal angle and reads 



n(a, (3) = + (n s — noo)\/2(cosha: — cos/3)x 

roo cosher cosh (2n - P)t . . . 

x / -i- -J—P_ 1/2 (cosh a) dr. (3.8) 

JO cosh tit cosh(vr - Q)t i /^+^ v 

Here n s is the density of the saturated vapor just above the liquid-air interface (or 
the potential of the conductor), noo is the ambient vapor density (or the value of the 
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potential at infinity), and P_i^2+ir( x ) are the Legendre functions of the first kind 
(despite the presence of it in the index, these functions are real valued). The surface 
of the lens is described by the two coordinate surfaces f3\ = n — 9 and = n + 9, 
and the f3 derivative is normal to the surface. Evaporation rate from the surface of 
the drop is therefore given by 

■/(«) = D±- d,n(a,P)\ 0=2n+pi = D ^_^l d,n(a, P)\^_ e , (3.9) 



/3=2tt+/3 1 R 



where D is the diffusion constant and hp = Rj (cosh a — cos f3) is the metric coefficient 
in coordinate f3. [Note that an incorrect expression for J with a plus sign in the metric 
coefficient was used in Eq. (A2) of Ref. (19).] Carrying out the differentiation, one 
can obtain an exact analytical expression for the absolute value of the evaporation 
rate from the surface of the drop as a function of the polar coordinate r: 



D(n s - rioo) 
J(r) = £ 



- sin 9 + \Pl (cosh a + cos 6>) 3 / 2 x 



roo cosher 

x / — tanh \(ir — 9)r] Pi /9 , • (cosh a) rdr 

JO coshvrr LV i/^+«tv 



(3.10) 



where the toroidal coordinate a and the polar coordinate r are uniquely related on 
the surface of the drop: 

R sinh a 

r = 1 n- I 3 ' 11 ) 

cosh a + cos 9 

Thus, expression (3.10) provides the exact analytical expression for the evaporation 
rate J as a function of distance r to the center of the drop for an arbitrary contact 
angle 9. This expression also corrects an earlier expression of Ref. (25) [Eq. (28)] 
where a factor of in the second term inside the square bracket is missing. 

The expression for the evaporation rate is not operable analytically in most cases, 
as it represents an integral of a non-trivial special function (which, in its turn, is 
an integral of some simpler elementary functions). In most cases, we will need to 
recourse to an asymptotic expansion for small contact angles 9 in order to obtain any 
meaningful analytical expressions. However, there is one exception to this general 
statement (not reported in the literature previously). An important quantity is the 
total rate of water mass loss by evaporation dM / dt, which sets the time scale for all 
the processes. This total rate can be expressed as an integral of the evaporation rate 
(defined as evaporative mass loss per unit surface area per unit time) over the surface 
of the drop: 

dM f „ , x r ——77 , . , f R 



dt 



-f J(r, <f>)yjl + (V/i) 2 rdrd(f) = - f J(r)yjl + (d r h) 2 2nrdr, (3.12) 
J ^4 J 



where the first integration is over the substrate area A occupied by the drop. This 
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expression actually involves triple integration: one in the expression above as an 
integral of J(r), another in expression for J(r) as an integral of the Legedre function 
of the first kind, and the third as an integral representation of the Legendre function 
in terms of elementary functions. However, one can significantly simplify the above 
expression and reduce the number of integrations from three to one. Substituting 

\Jl + (<9 r /i) 2 = 1/^1 — (r sin 9/R) 2 and evaporation rate (3.10) into Eq. (3.12) and 
changing variables and integration order a few times, one can obtain a substantially 
simpler result: 



— = -nRD(n s - noo, 
roo 1 + cosh 26t 



sin 9 
1 + cos 9 



A /-oo 1 + cosher 

+4 / tanh (vr - 9)r] dr 

JO smh27TT 



(3.13) 



[Note: Eq. (2.17.1.10) of Ref. (32) was used in this calculation.] This result together 
with the time derivative of expression (3.5) for total mass 

g£ = ir^fflg = * (3.14) 

dt dt (1 + cos#) 2 dt v 1 

provides a direct method for finding the time dependence of 9 for an arbitrary value 
of the contact angle. Combining the last two equations, one can obtain a single 
differential equation for 9 as a function of time t: 

d9 D(n s -noo) 2 [ sind 

— — " 77> -(I + cos 9) -+ 

dt pR 2 v ; [l + cos^ 

A roo 1 + cosh26'r , ,. nS . , ' 

+4 / tanh vr - 9)t] dr 

JO smh27TT 



(3.15) 



Having determined the dependence 9{t) from this equation, one can obtain the time 
dependence of any other quantities dependent on the contact angle, for instance, the 
time dependence of mass from relation (3.5), or any other geometrical quantities of 
the preceding section. 

In practice, however, we will always use the limit of small contact angles in all 
the analytical calculations of this and the subsequent chapters. This is the limit of 
most practical importance, and the usage of this limit will help us keep the physical 
essence of the problem clear from the unnecessary mathematical elaborations. While 
for numerical purposes our exact expressions are absolutely adequate, the analytical 
calculations in a closed form cannot be conducted any further for an arbitrary contact 
angle. 

Expanding the right-hand side of Eq. (3.15) in small 9, we immediately obtain that 
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the contact angle decreases linearly with time in the main order of this expansion: 

(3.16) 




where we introduced the total drying time tj defined in terms of the initial contact 
angle 6, = 9(0): 

tf = * PR% r- (3.17) 

In the main order, the total rate of water mass loss is constant and water mass also 
decreases linearly with time: 



« = ^l-f (318) 




This linear time dependence during the vast majority of the drying process was di- 
rectly confirmed in the experiments (21; 22); see Fig. 5.5. The dependence of the 
evaporation rate on radius (linearity in R) was also confirmed experimentally and is 
known to hold true for the case of the diffusion-limited evaporation (33). 

Before we entirely switch to the case of the small contact angles for the rest 
of this chapter, it is instructive to compare the small-angle analytical asymptotics 
of the preceding paragraph to the exact arbitrary-angle numerical results based on 
Eqs. (3.13) and (3.15). In Figs. 3.1 and 3.2, we plot the numerical solutions for 9(t) and 
M(t), respectively, for several values of the initial contact angle and compare them to 
the small- angle asymptotics of Eqs. (3.16) and (3.18). In these figures, Q{ and Mj are 
the initial contact angle and the initial mass of water in the drop, respectively. Note 
that in these graphs tf is not the total drying time for each Of, instead, it is just the 
combination of the problem parameters defined in Eq. (3.17), which coincides with 
the total drying time only when 9} — > 0. As is clear from these graphs, the total drying 
time increases with the increasing initial contact angle. However, it does not increase 
linearly as prescribed by the asymptotic expression (3.17); instead, it grows faster 
(Fig. 3.3). Fig. 3.3 demonstrates that the small-angle approximation works amazingly 
well up to the angles as large as 45 degrees (this can also be seen in Figs. 3.1 and 3.2). 
Therefore, working in the limit of small contact angles for the rest of this chapter, 
we will not loose any precision or generality for the typical experimental values of 
this parameter. Lastly, we note that the large-angle corrections may be responsible 
for the observed non-linearity of the experimentally measured dependence M(t), as 
is clear from the comparison of Fig. 3.2 (theory) and Fig. 5.5 (experiment). 

Expression for the evaporation rate (3.10) becomes particularly simple in the limit 
of small contact angles. Employing one of the integral representations of the Legendre 
function in terms of elementary functions [Eq. (7.4.7) of Ref. (23)], taking limit 9 — > 
and conducting a number of integrations, it is relatively straightforward to obtain the 
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Figure 3.1: Numerical results: Dependence of contact angle 9 on time t. Different 
curves correspond to different initial contact angles; values of parameter 9{ are shown 
at each curve. The analytical result [Eq. (3.16)] in limit 9{ — > is also provided (the 
solid curve). 




Figure 3.2: Numerical results: Dependence of water mass M on time t. Different 
curves correspond to different initial contact angles; values of parameter 9{ are shown 
at each curve. The analytical result [Eq. (3.18)] in limit 9^ — > is also provided (the 
solid curve). 
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Figure 3.3: Dependence of the total drying time on the initial contact angle. The 
solid curve demonstrates the analytical result for small contact angles [Eq. (3.17)]; 
the dashed curve represents the exact numerical result [inferred from Eq. (3.15)]. 



following result: 



D(n s -noo)2 a 
J{r) = 5 cosh- (0->O), (3.19) 



which, upon identification cosh a = (R 2 + r 2 )/(R 2 - r 2 ) for 6 = 0, can be further 
reduced to 

J « = ^7l=r <^ 0) - (3 - 20) 

Thus, for thin drops the expression for the evaporation rate reduces to an extremely 
simple result featuring the square-root divergence near the edge of the drop. The 
same result could have been obtained directly if we solved an equivalent electrostatic 
problem for an infinitely thin disk instead of the double-convex lens. It is particularly 
rewarding that after all the laborious calculations the asymptotic of our result is in 
exact agreement with elementary textbook's predictions [see Ref. (34) for the deriva- 
tion of the one-over-the-square-root divergence of the electric field near the edge of 
a conducting plane in the three-dimensional space]. Eq. (3.20) will be widely used 
below for all thin drops of circular geometry. 

For the sake of completeness, it is also interesting to note the opposite limit of the 
expression (3.10), when the surface of the drop is a hemisphere (0 = tt/2). In this 
limit, a similar calculation can be conducted [also employing Eq. (7.4.7) of Ref. (23) 
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and a couple of integrations], and the following simple result can be obtained: 

J(r) = D ^- n ^ (0 _> ^2). (3.21) 

R 

Thus, a uniform evaporation rate is recovered. Again, this result is in perfect agree- 
ment with the expectations; the same result could have been obtained if we directly 
solved the Laplace equation for a sphere (the hemispherical drop and its reflection in 
the substrate). The uniform evaporation rate is a result of the full spherical symmetry 
of such a system. Similar exact results can also be obtained for a few other discrete 
values of the contact angle (e.g. for 9 = 7r/4). 



3.3 Flow field 



With h and J in hand, we are in position to find the flow inside the drop according 
to prescription (2.12)-(2.13). In circular geometry this task is particularly easy, since, 
due to the symmetry, there is no component of the velocity (i.e. the flow is directed 
radially outwards) and v r is independent of <fi. Equations (2.12) and (2.13) can be 
combined to yield 



v r (r,t) = — l - f I - + d t h] rdr 
rhJQ \p ) 



(3.22) 



for the radial component of the velocity. Plugging Eqs. (3.3), (3.16), and (3.20) into 
the above equation, one can obtain the following expression for the flow inside the 
drop under assumption of small contact angle: 



v r (r,t) 



R 2 



( 



A(t f - t)r 



1- Ur 



(3.23) 



This is the final expression for the velocity that we were looking for. 

The velocity diverges near the edge of the drop with a one-over-the-square-root 
singularity in (R — r) at the contact line. This could have been deduced directly from 
the conservation of mass (2.3), where the divergent evaporation rate must be com- 
pensated by the divergent velocity (since the free-surface height is a regular function 
of coordinates and, moreover, vanishes near the contact line). Physically, change of 
volume near the edge becomes increasingly smaller as the contact line is approached 
and hence the outgoing vapor flux must be matched by an equally strong incoming 
flow of liquid. 

In addition, the velocity diverges at the end of the drying time. Since the amount 
of liquid removed by evaporation from the immediate vicinity of the contact line of 
a thin drop remains approximately constant, the amount of incoming fluid must also 
stay approximately constant over the drying time. This flux to the region in the 
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vicinity of the contact line is proportional to v r h, and hence velocity must scale as 
1/h in terms of its time dependence. Therefore, since the height decreases linearly 
with time (in the main order in 9), time dependence l/(tf — t) is to be expected for 
the velocity. The result above displays this expected behavior. 



3.4 Solute transfer and deposit growth 

Shape of the streamlines in the highly symmetric case of circular drops can be pre- 
dicted without any calculations: these are the straight lines from the center of the 
drop to the contact line along the radius. The parameter characterizing each of these 
streamlines is simply the polar angle <fi of the corresponding radial direction, and the 
variable along the streamline is r. The time it takes an element of fluid to reach 
the contact line can be inferred from the differential equation v r (r,t)dt = dr. Since 
dependences on distance and time in v r (r,t) of Eq. (3.23) are separable (factorize), 
the integration of one side of this differential equation from r\ (the initial location 
of an element of fluid) to R in distance and the integration of the other side from 
to t in time are straightforward to carry out. The resulting dependence of the time 
it takes an element of fluid to reach the contact line on the initial position of this 
element of fluid can be written as 



v) + 



R 



3/2 

= 1. (3.24) 



Clearly, t — when r\ = R, and t — tf when = 0, as expected from the physical 
intuition; however, the intermediate behavior is quite non-trivial, and no intuition 
would be of much help. 

The mass of the deposit accumulated at the contact line by time t is given by 
Eq. (2.15), which for the case of azimuthal symmetry can be rewritten as 

f R 

m = c h(r, 0) 2vrrrfr, (3.25) 
Jn 

where h(r, 0) refers to the expression for the equilibrium surface shape (3.3) at the 
beginning of the drying process (with 9 = 6i at t = 0). Simple integration yields 



m = rriQ 



1 - 



R 



(3.26) 



where tuq = ttcR^9i/A is the total mass of the solute present initially in the drop. 
When an element of fluid starts from a vicinity of the contact line (r^ R), the 
accumulated deposit mass is still virtually zero when it arrives at the contact line. 
When an element of fluid starts from the center of the drop (rj <^ R), virtually all 
the solute (m ~ uiq) is already at the contact line by the time the element of fluid 
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reaches it. 

Eliminating from expressions (3.24) and (3.26), one can finally obtain the de- 
pendence of the deposit accumulated at the contact line on drying time t: 



m = mo 



4/3 



(3.27) 



which can be rewritten in a more symmetric form 




1. 



(3.28) 



(We do not have a simple and intuitive explanation for this symmetry.) No solute 
(m = 0) is at the contact line at the beginning of the drying process (t — 0), and all 
the solute (m = mo) is accumulated at the contact line when the drop has completely 
dried (t — tf). At early times, the deposit mass scales with the drying time as a 
power law with exponent 4/3: 



/3t\ 4/3 

m^mo 



(t«t/). 



(3.29) 



As we will show in the next chapter, this scaling with time is universal for early drying 
times and should be observed in any geometry of the drop (as long as the contact 
line is locally smooth). At the end of the drying process, the rate of the deposit 
accumulation diverges as (tf — t) -1 / 4 : 



m ~ mo 




(\tf-t\<£t f ). 



(3.30) 



It is this final divergence that is responsible for the experimentally observed 100% 
transfer of the solute to the edge. 

Thus, in the case of circular geometry, it is relatively simple to obtain the desired 
scaling of the deposit mass at the contact line with time. The situation is far more 
complex in the case of angular geometry. This case is considered in the following 
chapter. 



CHAPTER 4 

DEPOSIT GROWTH FOR ZERO-VOLUME PARTICLES 
IN POINTED EVAPORATING DROPS 

In this chapter we will provide the asymptotic solution to the problem of solute 
transfer in the case of angular drops. In this geometry, the droplet is bounded by an 
angle a in the plane of the substrate (Fig. 4.1). We choose the origin of the cylindrical 
coordinates at the vertex of the angle, so that the angle occupied by the drop on the 
substrate is < r < oo and —a/2 < < a/2 (Fig. 4.2). Given the complexity of 
the geometry, we seek an approximate solution that captures the essential physical 
features and correct at least asymptotically. Here, in the geometry of an angular 
sector, the only possible locations of singularities and divergences (which normally 
govern properties of the solutions to differential equations) are at the vertex of the 
angle (i.e. at r = 0) and at its sides (i.e. at = ±a/2). Therefore, the most 
important physical features will be correctly reflected if asymptotic results as r — > 
and as <fi — > ±a/2 are found. Thus, we limit our task to determining analytically 
only the asymptotic power-law scaling for most quantities, and this task proves to be 
sufficiently challenging by itself. We will also provide some numerical results that do 
not rely on these assumptions. 

Most of the results presented in this chapter have been previously published in 
two our earlier papers (26; 27), and hence we will often refer to those works for further 
details. 



4.1 Surface shape 

The boundary problem for the equilibrium surface shape consists of the differential 
equation (2.11) and the boundary conditions at the vertex and at the sides of the 
angle: 

h(0, 0, t) = h(r, -a/2, t) = h(r, a/2, t) = 0. (4.1) 

Equation (2.11) represents the fact that the local mean curvature is spatially uniform, 
but changes with time as the right-hand side (Ap) changes during the drying process. 
The asymptotic solution to the boundary problem (2.11), (4.1) was found in our 
earlier paper (26). The result turned out to have two qualitatively different regimes 
in opening angle a (acute and obtuse angles) and can be written as 

Mr,M = g^. (4.2) 
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(b) 

Figure 4.1: (a) A water droplet with a sector-shaped boundary on the plane substrate 
(side view), (b) The same droplet pictured from another point (top view). Black lines 
are the grooves on the substrate necessary to "pin" the contact line. (Courtesy Itai 
Cohen.) 
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Figure 4.2: Geometry of the problem. The plane of the figure coincides with the 
substrate. 



Here Ro(t) = a/Ap and is the only function of time in this expression; exponent v 
has a discontinuous derivative at a = n/2 and is shown in Fig. 4.3; and 



1 / cos 20 _ i \ 
1 V cos a J 



if 0<a<§ (^ = 2), 

if \ < a < 7r (u — ix/a). 



(4.3) 



The constant C cannot be determined without imposing boundary conditions on h 
at some curve on the side of the drop furthest from the vertex of the angle. It is 
restricted by neither the equation nor the side boundary conditions, and thus, is not 
a universal feature of the solution near the vertex of the angle. The constant C can 
(and does) depend on the opening angle a. As we showed in the earlier paper, this 
constant must have the following diverging form near a = n/2: 



C 



1 



4a - 2% 



+ C + 0{a - tt/2) 



(4.4) 



where Cq is independent of a. We will adopt this form of C (with Cq set to unity) 
for all numerical estimates for obtuse opening angles. 

Two different values of v corresponding to the acute and obtuse angles give rise to 
the two qualitatively different regimes for surface shape. This difference can best be 
seen from the fact that the principal curvatures of the surface stay finite as r — > for 
acute angles and diverge as a power of r for obtuse angles. This qualitative difference 
can be observed in a simple experimental demonstration, which we provided in our 
earlier work (26). We refer to that earlier work for further details and discussion. We 
only note here that the asymptotic r — > at the vertex of the angle actually means 
r <^ Rq (which is typically of the order of a few millimeters for water under normal 
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Figure 4.3: Dependence of exponent v in the power law h(r) [Eq. (4.2)] on opening 
angle a, after Ref. (26). 



conditions), and that V/i is indeed small for r <^ Rq and can again be safely neglected 
with respect to unity (i.e. the free surface is nearly horizontal in the vicinity of the 
tip of the angle) as was asserted earlier. 

4.2 Evaporation rate 

As we explained in Chapter 2, in order to determine the evaporation rate one needs 
to solve an equivalent electrostatic problem for a conductor of the shape of the drop 
and its reflection in the plane of the substrate kept at constant potential. In geometry 
of the angular sector this shape resembles a dagger blade, and one has to tackle the 
problem of finding the electric field around the tip of a dagger blade at constant 
potential in infinite space. 

If one decides to account for the thickness of the blade [given by doubled h(r, 4>, t) 
of Eq. (4.2)] accurately, it becomes apparent that there is no hope for any analytical 
solution in this complex geometry. This is suggested by both the bulkiness of the 
expressions for a round drop with a non-zero contact angle in the preceding chapter 
and direct attempts to find the solution. However, taking into account that near 
the tip V/i is very small and hence the thickness of the blade itself is very small, 
we can approximate our thick dagger blade with a dagger blade of thickness zero 
and the same opening angle (i.e. with a flat angular sector). In the limit r — > the 
contact angle scales with r as (r/Ro) v ~ l and hence goes to zero. Thus, only the 
flat blade can be considered up to the main order in r. This approximation would 
not be adequate for determining the surface shape or the flow field, but it is perfectly 
adequate for finding the evaporation rate. We will discuss possible corrections to this 
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result later in this section. 

The problem of finding the electric field and the potential for an infinitely thin 
angular sector in the three-dimensional space requires introduction of the so-called 
conical coordinates (the orthogonal coordinates of the elliptic cone) and heavily in- 
volves various special functions. Luckily, it was studied extensively in the past (35; 
36; 37; 38), although the results cannot be expressed in a closed form. An impor- 
tant conclusion from these studies is that the r and dependences separate and that 
the electric field near the vertex of the sector scales with r as a power law with an 
exponent depending on the opening angle a: 



J(r, 0) oc r M_1 J(0). (4.5) 

Here 

(4.6) 



-L%(9*, 0*) = M/i + 1 W, 0*) (4.7) 



i cos 0* I 89* 

and fi and Y^{9*,(j)*) are the eigenvalue and the eigenfunction, respectively, of the 
eigenvalue problem 

with Dirichlet boundary conditions on the surface of an elliptic cone (degenerating 
to an angular sector as 9* — > 7r). In the last relation L 2 is the angular part of 
the Laplacian in conical coordinates (r, 9*, 0*). On the surface of the sector (i.e. at 
9* = 7r) the relation between the conical coordinate 0* and the usual polar coordinate 
is sin0 = sin(a/2) sin0*. We refer to work (37) for further details. Here we notice 
only that neither \x nor Y^O*, 0*) can be expressed in a closed analytic form; however, 
the exponent /i can be computed numerically and is shown in Fig. 4.4 as a function of 
a. Note that this exponent is lower than similar exponents for corresponding angles 
for both a wedge (a two-dimensional corner with an infinite third dimension) and a 
circular cone, which are also shown in the figure. Both these shapes (wedge and cone) 
allow simple analytical solutions (34) but none of them would be appropriate for the 
zero-thickness sector. 

Despite the unavailability of an explicit analytical expression for J(0), its analytic 
properties at = and at = ±a/2 are quite straightforward to infer. Indeed, J(4>) 
is an even function of 0; therefore, J'(0) = (as well as any other odd derivative on 
the bisector) and 

J(0) = J(O)+O(0 2 ) (4.8) 

for small 0. Obviously, J(0) is positive. On the other hand, at = ±a/2 the leading 
asymptotic of the evaporation rate (or the electric field) is known to be (A0) -1 / 2 
with exponent —1/2 corresponding to the edge of an infinitely thin half plane in the 
three-dimensional space (34). (We have introduced notation A0 = a/2 — |0| in the 
previous line.) If one were to correct this asymptotic in order to reflect the non-zero 
contact angle at the edge of the sector, the asymptotic form at = ±a/2 would have 
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Figure 4.4: Solid line: Dependence of exponent \x in the power law J(r) [Eqs. (4.5) 
and (4.12)] on opening angle a for an angular sector, after Refs. (36; 37; 38). The 
same dependencies for a wedge (dashed line) and a cone (dotted line) are also shown. 



to be written as 

J(0) » J*(A0)- A 
where J* is a positive constant and 

71-29 



A 



2u-29' 



(4.9) 



(4.10) 



This result corresponds to the divergence of the electric field along the edge of a 
wedge of opening angle 29 [both the drop and its reflection contribute to the opening 
angle, hence a factor of 2] (34). However, accounting for the non-zero 9 is a first-order 
correction to the main-order result A = 1/2. This can be seen from the expression 
for the contact angle: 



9 = arctan 



v-l 



~ti {a/2) 



v-l 



(X 



(4.11) 



For all opening angles v > 1 (except a = rr where v = 1). Thus, the correction due to 
the non-zero contact angle can indeed be neglected in the main-order results, and A 
should indeed be set to 1/2. Nevertheless, we will keep the generic notation A for this 
exponent in order to keep track of the origin of different parts of the final result and 
in order to account properly for the case a = it in addition to the range of opening 
angles below n. The numerical value of A will be assumed to be 1/2 in all estimates. 
The analytical results below will employ the asymptotics of this paragraph. 
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Thus, we will use the following expression for the evaporation rate J: 

J(r ' 0) = J °(^) M J(0) ' (412) 

where function J(4>) is defined in Eq. (4.6) with asymptotics (4.8) and (4.9). Here 
we broke down the constant prefactor into two pieces: a distance scale \fA (where 
A is the substrate area occupied by the drop) and all the rest Jq (which is of the 
dimensionality of the evaporation rate). Trivially Jq is directly proportional to the 
product of the diffusion constant and the difference of the saturated and the ambient 
vapor densities D(n s — noo), as was the case in the round-drop geometry. 

The evaporation rate above does not depend on time and the same form of J 
applies during the entire drying process, since the diffusion process is steady and the 
variation of the contact angle with time does not influence the main order result (4.12). 
The same is true for the total rate of mass loss dM/dt since 



dM 
~df 



- I J\[\ + (Vh) 2 rdrd(f) « - f J rdrdcj) oc -J A, (4.13) 



where the integrations are over the substrate area occupied by the drop. We saw 
the same behavior in circular geometry. The constancy of this rate during most of 
the drying process was also confirmed experimentally (21). This fact can be used 
to determine the time dependence of the length scale Rq of Eq. (4.2) [and hence of 
the pressure pq of Eq. (2.11)] explicitly, as the mass M of a sufficiently thin drop is 
inversely proportional to the mean radius of curvature Rq: 

nA 2 

Moc^-, (4.14) 
Kq 

where we retained only the dimensional quantities and suppressed all the numerical 
prefactors sensitive to the details of the drop shape. From the last two equations one 
can conclude that 

x(l^l«-4 (4-15) 
dt\Ro) pA v ; 

and remains constant during most of the drying process. Hence, 

*bM - J^-, (4.16) 

where Rqi is the initial mean radius of curvature (Rqi = Rq(0)) and tf is the total 
drying time: 

t f oc (4.17) 
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Thus, at early drying stages (t <C tf) scale Rq grows linearly with time; this time 
dependence will be implicitly present in the results below. However, it is very weak at 
sufficiently early times and will be occasionally ignored (by setting Rq Rqi) when 
only the main-order results are of interest. 



4.3 Flow field 



As in the case of circular geometry, with h and J in hand, we proceed by solving 
Eq. (2.12) for the reduced pressure tp. Assuming power-law divergence of ip as r — > 
and leaving only the main asymptotic (which effectively means that we neglect the 
regular term dth with respect to the divergent one J/p), we arrive at the following 
asymptotically-correct expression for ip: 

1>(r, 0, t) = - —3i (4-18) 

where time-dependence is implicitly present via Rq and the function ip(<f>) is a solution 
to the following differential equation: 

/l(0) ^(0) 

(the combination (// + l)(3z/ — // — 1) is positive for all opening angles). Computing 
v according to prescription (2.13), we obtain the depth-averaged flow field 

v = Vr f + v^> (4.20) 

with components 



v r (r, 0, t) = -(31/ -iJL- 1)^——^ /^(0)</,(0) (4.21) 

P y/Af* Rq V+1 

and 

v<p(r, 0, t) = ^—^-—^(0)^(0). (4.22) 

Thus, one needs to solve Eq. (4.19) with respect to -0(0) in order to know the flow 
velocity. 

We were not able to find an exact analytical solution to this equation; however, 
we succeeded in finding approximate solutions on the bisector (|0| <C a/2) and near 
the contact line (A0 <C a/2), which represent the two opposite limits of the range 
of 0. (Again, we define A0 = a/2 — |0|.) We describe these two solutions in great 
detail in our earlier work (27), and since the solutions themselves are not needed in 



33 



this abbreviated account, we refer to this earlier work for further details. Here we 
will only mention the asymptotics of both solutions (which can actually be inferred 
directly from the differential equation, without even solving it). Near the contact 
line, in limit A<p — > 0, the asymptotic of the solution is 

^(0) oc J * (A0) ~ A ~ 1 (4.23) 
(1-A 2 )|// (a/2) | 

On the bisector, in the opposite limit 0^0, the asymptotic is 

\ (4.24) 



4>(<f>) oc 4>(o) + x - 



The reduced pressure on the bisector -0(0) is positive. The value of -0(0) cannot be 
determined from the original differential equation; one needs to employ an integral 
condition resulting from the equality of the total influx into a sector of radius r by 
flow from the outer regions of the drop and the total outflux from this sector by 
evaporation: 

ra/2 rr ra/2 

p J \v r \hrd^) — J ^Jrdrdcj). (4-25) 

This condition is similar to an analogous condition for the circular drop (which was 
written for the entire drop instead of only a certain part of it, since, unlike the infinite 
sector, the round drop was finite). Upon simplification this condition reduces to the 
following equation defining the constant prefactor -0(0): 



"/ 2 r . ,.r 3 , 



[(pi + 1)^(0)^(0) - J(0)J dcf) = 0. (4.26) 

Obviously, -0(0) is proportional to J(0)/i~ 3 (0). 

In order to compensate for the unavailability of the exact analytical solution to 
Eq. (4.19) we also approached this problem numerically. We will not describe the 
numerical procedure in full here as it was described in great detail in our earlier 
paper (27). Here, we will only mention that we use two different trial forms of J(4>) 
that are simple and at the same time satisfy proper asymptotics (4.8) and (4.9). 
These model forms are neither exact nor the only ones satisfying asymptotics, but 
they allow one to avoid the numerical solution of the eigenvalue problem (4.7) and 
thus not to repeat the elaborate treatment of works (35; 36; 37; 38). We use the 
results based on these trial forms only for illustrative purposes in order to picture 
general behavior of the solution for arbitrary values of the argument (not only in the 
limiting cases). As we discussed in work (27), the difference between the numerical 
results based on the two model forms of J(<j>) did not exceed 10-15% in most cases, 
and we have all reasons to believe that at least the orders of magnitude obtained by 
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Figure 4.5: Typical behavior of the numerical solution and the analytical asymptotics 
of function if>((f)) for two values of opening angle. 



this approximation are correct. We would like to emphasize that only the numerical 
graphs based on the choice of J((f>) are affected by these simplified forms; all the 
analytical results below do not rely on a particular form of J((f>) and use only the 
analytical asymptotics of the preceding section. 

The numerical solution to Eq. (4.19) satisfying conditions (4.26) for if)(0) and 
■0(0) = for ip'(0) was found for the two model forms of J ((f)) and for approximately 
20 different values of the opening angle. In all cases perfect agreement between 
the numerical solution and the analytical asymptotics of this section was observed. 
Two examples of the numerical solution together with the analytical asymptotics are 
provided in Fig. 4.5 for opening angles 70° and 110°. 

Characteristic behavior of the velocity field (4.20) is shown in Fig. 4.6 for a = 70° 
and a = 110°. Note that despite the fact that the exponent (/i — v) of the power law 
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in r is not a smooth function of a, the qualitative behavior of the flow field does not 
visibly change as the opening angle increases past the right angle. 

The velocity diverges near the edge of the drop. At the sides of the angle, only 
the component diverges, as can be seen from Fig. 4.6 and expression (4.22). This 
divergence near the sides is of exactly the same one-over-the-square-root dependence 
on distance as the divergence in the circular drop and has exactly the same origin. 
In addition, there is a divergence of both components of the velocity as r^~ v when 
r — > 0, as is apparent from expressions (4.21) and (4.22). This divergence is entirely 
new and is due to the presence of the vertex. The exponent of this divergence depends 
on the opening angle as both fi and v do. This is the first in a set of indices for the 
pointed drop that are universal but depend on the geometry. Similar indices will also 
be encountered for other physical quantities throughout this chapter. 



4.4 Streamlines 



r(0) = rg exp 



is, -„-!)/; 



(4.27) 



One must know the shape of the streamlines in order to be able to predict the scaling 
of the deposit growth, and now, with velocity field in hand, we have everything needed 
to compute it. Integration of the velocity field (4.21)-(4.22) according to Eq. (2.14) 
yields the streamline equation, i.e. the trajectory of each particle as it moves with 
the fluid: 

ra/2 j(g)rfg - 
4>'(0 J 

where we assume that is positive here and everywhere below (the generalization to 
the case of negative is obvious as all functions of are even). Thus, r = rg when 
= a/2, so that rg is the distance from the terminal endpoint of the trajectory to 
the vertex. In limit A0 — > the integral in the exponent goes to zero (quadratically 
in A0) and the streamline equation reduces to 

r ^ r (4.28) 

[see Ref. (27) for details of all calculations in this section]. The streamlines are 
perpendicular to the contact line (up to the quadratic terms in A0). This is in 
good agreement with what one would expect near the edge of the drop, since the 
azimuthal component of the fluid velocity diverges at the side contact line while the 
radial component goes to zero. In limit — » the integral in the exponent diverges 
as ln[(a/2)/0] and hence the result reads 

(4.29) 
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Figure 4.6: Characteristic behavior of flow field for two values of opening angle. Each 
arrow represents the absolute value and the direction of velocity v at the point of 
arrow origin. 
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Here we introduced exponent e: 



1 



(4.30) 



e = 



(// + !)(!-/)' 



where 



a/2 (h^)\ S j(<t>) 



U(q)7 m 



d(f) 




(4.31) 



As we showed in Ref. (27), < / < 1 for all a. The positiveness of exponent e 
follows both from this fact and from the fact that the trajectory r(0) necessarily has 
to diverge as — > (as solute comes from the outer regions of the drop). 

Note that the last expression for parameter / is independent of the prefactors 
of each function of 0, particularly, it is independent of the evaporation intensity 
Jq and the constant non-universal prefactor C of the surface shape, reinforcing the 
general conclusion of Chapter 2 that the shape of the streamlines is universal for each 
geometry. The same conclusion can also be reached directly from Eq. (4.27), which 
is independent of Jq and C and thus is universally correct regardless of the overall 
intensity of evaporation and the non-universal features of the drop boundary. 

We cannot compute I and e explicitly, since we do not have an analytical expres- 
sion for ■0(0). However, we can gain some idea of the behavior of these parameters by 
using approximate forms of J{4>) as we did in the preceding section. Fig. 4.7 demon- 
strates the characteristic behavior of exponent e as a function of the opening angle, 
obtained numerically on the basis of the two model forms for function J ((f)). In order 
to obtain this plot, equation (4.19) was solved with respect to i/j(4>) numerically for 
each a, and then expression (4.31) for I was employed. As can be observed in this 
graph, the two model forms of J(<f>) lead to the plots of very similar shape, but shifted 
by no more than 10% in the whole range of the opening angles. Thus, we conclude 
that Fig. 4.7 provides a good estimate for the qualitative behavior and the order of 
magnitude of exponent e(a). 

Typical shape of the streamlines is shown in Fig. 4.8 for a = 70° and a = 110°. It 
was based on one of the model forms for function J (</>), and involved the corresponding 
numerical solutions for function -0(0) (Fig. 4.5) employed in Eq. (4.27). This shape 
is practically insensitive to the model form of J(0), and almost an identical copy of 
this graph was obtained for the other model form. 

The distance from a point on a streamline to the bisector scales with as 0r(0) oc 
1_e when 0^0. Since e > 1 (Fig. 4.7), this distance increases when decreases. 
Thus, the streamlines diverge away from the bisector when 0^0, and hence they do 
not originate on the bisector. An incoming element of fluid initially located close to 
the bisector moves towards this bisector, reaches a minimum distance and then veers 
away towards the contact line. One can also arrive at the same conclusion having 
started from the ratio of velocity components. As can be shown [see Ref. (27)], for 
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Figure 4.7: Dependence of exponent e in the power law r((f>) [Eq. (4.29)] on opening 
angle a. The two curves correspond to the two model forms for function J(0). 



small (in limit — > 0) the ratio of the velocity components is tu/tv = — 0/e- This 
ratio represents the angle between a streamline and a coordinate line = 0o at any 
point (0o 5 r (0o)) on that streamline. Since e > 1, the absolute value of this angle 
is less than \(f>\, and therefore, despite the opposite sign of this angle, the streamline 
diverges away from the bisector for small 0. This tendency can also be observed 
directly in Fig. 4.8. 

Another feature apparent from Fig. 4.8 is the self-similarity of all the streamlines. 
As is clear from equation (4.27), the only scaling parameter of the family of streamlines 
is ro, and therefore all the streamlines can be obtained from a single streamline (say, 
the one with tq = 1) by multiplying its r-coordinate by different values of tq. 

4.5 Solute transfer and deposit growth 

The time it takes an element of fluid to reach the contact line can be inferred from 
either v r (r, 0, t)dt = dr or f^(r, 0, t)dt = r d(p on a streamline (i.e. with relation (4.27) 
between r and 0). The only term with time dependence in expressions (4.21) and 
(4.22) is i?Q~ . As we mentioned earlier, the time dependence of this term is weak 
and can be safely ignored for times t <C tt, where tt is the total drying time defined 
in Eq. (4.17). It will be seen later in this section that these times are the only times 
of interest in the geometry of the pointed drop since everything that happens closer 
to the end of the drying process is not universal and depends on the details of the 
boundary shape outside the sector of interest (i.e. on the shape of the contact line 
on the side of the drop furthest from the vertex). Therefore, restricting our attention 
to only the universal features of the solution, we can consider only the times t <C 
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Figure 4.8: Typical shape of the streamlines for two values of opening angle (for the 
same values as in Figs. 4.5 and 4.6). 
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and ignore the dependence of Rq (and hence v r and va) on time. Then the time and 
the coordinate dependences in the differential equations above separate trivially and 
the time it takes an element of fluid to reach the contact line can be calculated as 

W2 rdd) fro dr 
t= —= \ — , (4.32) 

Hi Vff, Jn v r 

where (r^, 0j) are the initial coordinates of that element of fluid. Taking into account 
expressions (4.21) and (4.22) for the velocity components and relation (4.27) between 
r and on a streamline, one can reduce both integrals to 



t=t °L 



«/2 ex P 



;,-, + D(3,-,-i)/;W 



#, (4.33) 



h\ct>W{ct>) 

where to is a combination of system parameters with the dimensionality of time: 

to = j-VA! i - 1 ltf+ 1 r v -'* 1 (4.34) 

and Rqi = -Ro(O). Within this time all the solute that lays on the way of this element 
of fluid as it moves toward the contact line becomes part of the deposit (highlighted 
area in Fig. 4.9). The mass dm of this deposit (accumulated on the contact line 
between tq and tq + dro) can be found by integrating h(r, <p, 0) of Eq. (4.2) over 
area dA swept by this infinitesimal volume and multiplying the result by the initial 
concentration c of the solute [see Eq. (2.15)]. Employing relation (4.27) once again, 
we obtain: 



dm = c r ^^ r /2 M0)expL + 2X3^-1) J 



R 0i 



*/2 gQ< 

4>'(0 . 



d(j>. (4.35) 



Dependence dm(t) can now be found by eliminating 0j from results (4.33) and (4.35). 

Exact analytical calculation of the dependence m(t) is not possible for an arbi- 
trary starting point (r^/pj) on a streamline since no analytical expression for "0(0) 
is available for arbitrary and since integrals in Eqs. (4.33) and (4.35) cannot be 
computed analytically for arbitrary even if ip(<p) were known. However, there are 
two important cases that can be tackled analytically: (a) early times, when the initial 
point is close to the contact line (i.e. when A<^ < a/2 and pa tq) and only the 
solute between that initial point and the contact line is swept into the edge deposit 
(the starting point is on isochrone 1 of Fig. 4.9 or closer to the contact line), and (b) 
intermediate times, when the initial point is close to the bisector (i.e. when |0^| <^ a/2 
and ri ^> vq) and virtually all the solute between the bisector and the contact line is 
swept into the edge deposit (the starting point is on isochrone 4 of Fig. 4.9 or further 
from the vertex). Situations between these two limiting cases (highlighted area in 
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Figure 4.9: Qualitative sketch: mutual location of streamlines (the two lines with 
arrows) and isochrones (the four numbered lines). Solute moves along the streamlines 
towards the contact line (the bold line). Shaded area is swept by an infinitesimal 
element of fluid between the two infinitesimally close streamlines as that element 
moves towards the contact line. The isochrones are the geometric locations, starting 
from which the solute reaches the contact line at the same time. Solute from isochrone 
1 reaches the contact line first; solute from isochrone 4 reaches the contact line last. 



Fig. 4.9 demonstrates one of them, starting points on isochrones 2 and 3 would cor- 
respond to some other) can be extrapolated on the basis of continuity of the results. 
Since our region is indefinitely smaller than the drop as a whole, we may treat regions 
(a) and (b) assuming that a negligible fraction of the drop has evaporated. At some 
later stage that we call the late-time regime, an appreciable fraction of the drop has 
evaporated, and the fluid trajectories have reached back into the bulk of the drop. 
In this late regime our asymptotic treatments are clearly not adequate to describe 
the flow as we did not specify any details of the drop geometry far from the vertex. 
Thus we cannot treat this regime by our methods, and only the properties of drying 
process at early and intermediate stages can be found from information in hand. 

Apart from the definitions based on the trajectories, the three regimes can be 
equivalently defined in terms of time t: (a) early times: t <C to; 3 ) intermediate 
times: to ^ t <^ t r, and (c) late times: t ~ tf. Here to is the characteristic time 
defined in Eq. (4.34) (this characteristic time depends on tq) and tf is the total drying 
time defined in Eq. (4.17). The equivalence of the definitions in terms of the initial 
position on a trajectory and in terms of time can be seen from equation (4.33). As we 
explained in the preceding paragraph, only the early and the intermediate times are 
of our interest as being independent on the rest of the boundary shape, and therefore 
the assumption t < employed at the beginning of this section is coherent with our 
intention to determine only the universal features of the results. 

As is clear from the definition of the intermediate-time regime, the necessary condi- 



42 



tion for its existence is tg ^ ~tf, which can be reduced to (rg/ v^4) 3 ^( r o/ Roi) 1 ' ~ 2 <S 1 
by combining Eqs. (4.34) and (4.17). For acute opening angles the latter condition is 
always true (27), and hence the intermediate regime is well defined for acute open- 
ing angles. For obtuse opening angles the condition above can be transformed into 
(ro/v^4)^~' i+1 <C 6j~ l where Q{ is the typical initial value of the contact angle in 
the bulk of the drop, i.e. far away from the vertex. This condition should be expected 
to be satisfied for not too small initial values of the bulk contact angle Q{. The closer 
to the vertex the trajectory endpoint is, the better this condition is obeyed. On 
the other hand, this condition is obeyed worse for larger opening angles. At exactly 
a = it the intermediate-time regime is indistinguishable from the early-time regime, 
and hence should not exist (27). 

At early times, the deposit growth is entirely due to the transfer of particles 
originally located near the contact line. Evaluating the integrals in expressions (4.33) 
and (4.35) in limit A0j — > 0, expressing A0j in terms of time and then substituting 
the result into the expression for mass, we obtain the mass of the deposit as a function 
of time: 



dm,. , r v Q +l 



h! (a/2) 



1 + A J* t 



2 



— (t >n ,)« c " ' - , — - . (4.36) 



dr " u/ 2 \1-A 



h' (a/2) 



to 



Note that to also depends on tq. Thus, at early times the deposit grows in time as a 
power law 

^(t,n.)«* J /< 1+ ^, (4.37) 

where the Tq arises from the Tq +1 prefactor and from the rg-dependence of to- Using 
Eq. (4.34), we find 

P = _(l-A)(l + -)-2/. (4 3g) 

1 + A 

and plot it in Fig. 4.11 as a function of opening angle (the early-time curve). Note 
that this result is independent of function ip(<f)) and hence precise as long as the factor 
J* in the asymptotic of the evaporation rate (4.9) is known. 

There are two important conclusions to be drawn from this result. One is that 
the power-law exponent of time 2/(1 + A) = 4/3 is exactly the same as in the case 
of a round drop considered in the preceding chapter. This should be of no surprise 
since close to the side of the angle (as well as close to the circumference of a round 
drop) the contact line looks locally like a straight line, and the solute "does not know" 
about the vertex of the angle or the curvature of the circumference. This exponent is 
determined entirely by the local properties of an infinitesimal segment of the contact 
line of length drQ and is independent of larger geometrical features of the system. 

The value 2/(1 + A) = 4/3 of the exponent of time can be obtained from a 
very simple argument, relying only on the assumptions that (a) the contact line 
is straight, (b) the streamlines are perpendicular to the contact line, and (c) the 
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Figure 4.10: An illustration of the derivation of the four-thirds law at a straight 
segment of the contact line. The contact line is normal to the plane of the figure. 
Length drQ is along the contact line and hence not shown. The flow is in the plane 
of the figure from left to right. 



distribution of the solute is uniform. Indeed, the mass of both the water and the 
solute is proportional to the volume of an element of fluid near the contact line 
(Fig. 4.10): dm oc (Al) 2 dr . All this mass should be evaporated from the surface 
of this volume element in some time t. The evaporation rate (per unit area) scales 
as J oc (AZ)~ A and therefore the full rate of mass loss is JdA oc (Al)~ x+l drQ. The 
time it takes this volume to evaporate can now be found as the ratio of its mass to 
the rate of mass loss: t = dm/(JdA) oc (AZ) 1+A Thus, (A/) oc tV(l+A) and hence 
dm/dr oc (A/) 2 oc t 2 /( 1+A ) as asserted. 

The other observation is the dependence on tq. First of all, the dependence 
on distance to the vertex itself is entirely new compared to the round-drop case 
(for obvious reasons). Second, since exponent (5 is always between —1 and 0, the 
singularity in tq is always integrable at tq = 0. Physically, this corresponds to the 
statement that the vertex of the sector does not dominate the sides and that the 
deposit accumulation at the vertex is not qualitatively different from the deposit 
accumulation on the sides. 

Similar calculation can also be conducted analytically for the intermediate times 
as well. The starting point of a streamline in the intermediate-time regime lies near 
the bisector, and the appropriate limit is 0^ — > 0. Calculating the time it takes an 
element of fluid to reach the contact line [Eq. (4.33)] and the mass accumulated at 
the contact line between tq and tq + dr-Q for this time [Eq. (4.35)] in limit fa — > 0, 
and then eliminating fa from the two results, we arrive at the dependence of mass on 
time: 

^ (t ' r ° ) ~ £ Xri (,-^ + l) + (^ + l)/ M0) 2 X 
/ t \ i+0*+i)V("-m+i) 

x Uu-n+ l)(Zv-n- l)h 2 (0)ip(0) — j . (4.39) 

Taking into account that a l so depends on rg, we finally conclude that the deposit 
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mass grows as a power law 



^(t,r )cxt 5 rl (4.40) 



where we introduced notations for the exponent of time 

s=1+ (t±jv (441) 

(v-n + 1) 

and for the exponent of tq, originating from both the prefactor Tq +1 and the rg- 
dependence of tQ [Eq. (4.34)]: 

7 = /i-(/i + l)/. (4.42) 

Unlike the early-time limit, this result depends on the entire functions J ((f)) and ip(4>), 
not only their asymptotics, via parameter I defined in Eq. (4.31). This is quite natural 
because an element of fluid passes through the entire drop when it moves from the 
bisector to the contact line, and hence knowledge of all quantities in the intermediate 
points is required. In the early regime, only the asymptotic properties near the edge 
are of importance. 

The exponent of time stays greater than one in the intermediate-time regime. 
Thus, the rate of mass accumulation dm/ dt continues to grow with time in this regime, 
and the deposit mass grows faster and faster. This result has a simple explanation 
for both the early- and the intermediate-time regimes. Since the initial distribution 
of the solute is uniform, and since the solvent evaporates, the solute concentration at 
any given volume increases with time. Thus, even though the fluid and the particles 
move along the same streamlines in practically constant velocity field (assuming that 
Ro(t) ~ Rqi at sufficiently early stages), the rate of mass accumulation also increases 
with time, since portions of solution arriving at the contact line at approximately 
constant rate have higher and higher solute concentration. Note that this mechanism 
and this general conclusion are in good agreement with the exact analytical result 
for the circular geometry that the rate of mass accumulation must diverge at the end 
of the drying process (as t — > tj) and that all the deposit must accumulate at the 
contact line by t = tj. 

Another observation is related to the exponent of tq. Since I < 1, then 7 > — 1. 
Therefore, the mass is integrable at rg = 0, and the statement of the early-time 
regime that the deposit accumulation at the vertex is not qualitatively different from 
the deposit accumulation on the sides continues to hold in the intermediate-time 
regime as well. Trivially, 7 < /i, as / > 0. 

The exponent of rg must be identically zero at any time for the opening angle of 
exactly a = ir. Indeed, at a = n the contact line is just a straight line (i.e. there 
is no angle at all), and therefore there is a full translational symmetry with respect 
to which point of this line should be called "vertex." Thus, the choice of rg = is 
absolutely arbitrary, and there can be no dependence on rg whatsoever. 
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Figure 4.11: Exponent of distance rg in the power law dm/dro(rQ) [Eqs. (4.37) and 
(4.40)] as a function of the opening angle for the two time regimes. The early-time 
curve corresponds to the exponent (5 of Eq. (4.38); the intermediate-time curves 
correspond to the exponent 7 of Eq. (4.42). The two curves for the intermediate-time 
exponent correspond to the two model forms for function J(<f>). 



Indices 7 and 5 are plotted in Figs. 4.11 and 4.12, respectively, as functions of the 
opening angle (the intermediate-time curves). The two intermediate-time curves on 
each graph correspond to the same two model forms for function J(<f>) as we used in 
Fig. 4.7. To facilitate the comparison of the results, we plot the exponents for the 
early- and the intermediate-time regimes in Figs. 4.11 and 4.12 together. 

The intersection of the exponents near a = n on both graphs can be attributed 
to a couple of reasons. First, it should be kept in mind that the plotted results for 
7 and 5 are based on a relatively arbitrary choice of model forms for function J(<f>), 
which, as we suspect, become increasingly inaccurate for large a. Second, as was 
explained just after Eq. (4.11), at exactly a = n the contact angle 6 is not small 
even for r <C Rq, and the correction to the exponent A due to this contact angle [see 
Eq. (4.10)] is comparable to the value 1/2 assumed in all numerical estimates. All in 
all, we believe that this intersection of the early- and intermediate-time exponents is 
an artifact of our formalism and should not be observed in reality, since the results 
for the two time regimes should be identical at exactly a = ir. At exactly a = it the 
exponent of rg must be equal to zero at any time and the exponent of time must be 
equal to 2/(1 + A) =4/3 at any time. 

In addition to the early- and the intermediate-time analytical asymptotics above, 
we also find the numerical solution for d?m/ 'dtdro(t). We find the time derivative of 
dm/dr-Q instead of dm/dr-Q itself in order to demonstrate the amount of mass arriving 
at the contact line at time t rather than the total mass accumulated by the time t. 
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Figure 4.12: Exponent of time t in the power law dm/dro(t) [Eqs. (4.37) and (4.40)] 
as a function of the opening angle for the two time regimes. The early-time curve 
corresponds to the exponent 2/(l + A) = 4/3; the intermediate-time curves correspond 
to the exponent 5 of Eq. (4.41). The two curves for the intermediate-time exponent 
correspond to the two model forms for function J(4>). 



We employ the chain rule to obtain d^m/dtdr-Q on the basis of Eq. (4.33) for t(<pi) 
and Eq. (4.35) for dm/ dro(4>j): 

d_ ( dm \ = 4j fe) = 
dt{dr )- *. 



c r v+ \ r 

to Bfc U 



ra/2 fafidf 

^(0 J 



(4.43) 



then use the numerical result for ip{4>) (Fig. 4.5) in order to find t(</>j) [Eq. (4.33)] 
and d*m/ dtdr§((f)j) [Eq. (4.43)] numerically, and finally create a log-log parametric 
plot d 2 m/dtdro vs. t, as shown in Fig. 4.13. The two curves in Fig. 4.13 correspond 
to the two values of a we used earlier (70° and 110°). As all numerical results of this 
section are, the plot is based on one of the model forms for function J ((f)), but very 
insensitive to the particular form of this function. This plot clearly demonstrates two 
different slopes (and hence two different time regimes) of each curve. The crossover 
between the two regimes (slopes) occurs around time t ~ to (^- e - near l n (^Ao) ~ 0)> 
and the early-time slopes are equal for both values of the opening angle (and equal to 
2/(1 + A) — 1 = 1/3 as to be expected from our early time results). All these numerical 
results are in excellent agreement with our analytical predictions, and the numerical 
values of time exponents compare very well with those of Fig. 4.12 (which should be 



47 



2-< 



70' 




-2 



-6 



-4 



-2 







2 



4 



6 



ln(t/t () ) 



Figure 4.13: Log-log plot of the numerical solution d m/dtdro(t) for two values of 
opening angle (for the same values as in Figs. 4.5, 4.6 and 4.8). 



corrected by —1 due to the differentiation with respect to time in Fig. 4.13). 

In a similar fashion we obtain a log-log plot for d^m/dtdrQ as a function of tq. 
We fix t, then express tq in terms of <pi by combining Eqs. (4.33) and (4.34), and 
finally determine ro(<^) and d 2 m/dtdro(<f)j) [Eq. (4.43)] numerically on the basis 
of the numerical result for ip{<p) of Fig. 4.5. The resulting log-log parametric plot 
d?m/dtdro(ro) is shown in Fig. 4.14 for the two values of the opening angle (70° 
and 110°). The purpose of this graph is to provide a snapshot of the deposit growth 
at any given moment of time t. For small tq the accumulation of the solute at the 
contact line is already in the intermediate-time regime, while for large tq the growth 
is still in the early-time regime. The threshold between the two regimes is defined by 
t = tQ. This condition can be reversed by solving Eq. (4.34) with respect to rg. The 
resulting value 



defines the threshold in terms of tq (at any moment of time t): the early regime 
corresponds to rg r*, and the intermediate regime corresponds to ro C r*. As 
can be seen from the numeric plot, the regimes indeed switch at tq r* {i.e. near 
ln(ro/r*) m 0). The intermediate-time slopes are almost equal for both graphs since 
the intermediate-time exponent 7 (the upper curves in Fig. 4.11) varies very weakly 
with a (7 ~ —0.135 for a = 70°, and 7 ~ —0.111 for a = 110°). Again, the numerical 
results are in excellent agreement with the analytical asymptotics, and the numerical 
values of exponents compare very well with those of Fig. 4.11. 

Figure 4.14 suggests that the rate of increase of dm/drQ has a sharp change of 
behavior as a function of tq for any given time. For small tq (intermediate times) 
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Figure 4.14: Log-log plot of the numerical solution d 2 m/dtdro(rQ) for two values of 
opening angle (for the same values as in Figs. 4.5, 4.6, 4.8 and 4.13). Parameter r* 
is defined by Eq. (4.44). 



function d mj ' &tdr§{r§) varies weakly, while for large rg (early times) it falls off more 
dramatically with increasing tq. The crossover point r* [Eq. (4.44)] moves outwards 
as a power of time t with exponent — /i + 1). (Note that this exponent involves 
only the accurately known functions of a.) This crossover point and its outward 
motion provide a clear-cut signature of our mechanism, and this signature should 
be the strongest for small opening angles (as Fig. 4.11 suggests). (In order to avoid 
possible non-universal effects from late times, one needs to measure the system before 
the late-time regime.) 

Probably, the most exciting feature of the angular-sector solution is its dependence 
on the opening angle. Unlike the round-drop case, there is an extra free parameter 
of the problem — the opening angle of the sector. All the results, including the 
exponents of the power laws, depend on this opening angle. Note that these exponents 
are universal, i. e. they do not depend on any other parameters of the system, except 
for the opening angle. At the same time, the only parameter they depend on is 
extremely easy to control — preparing an evaporating drop one can adjust the opening 
angle of the contact line at his will without any technical elaborations. Thus, for 
example, by suitably choosing the opening angle (and the time regime), one can 
create a predetermined power-law distribution of the solute along the contact line 
with virtually any exponent of distance between —1 and (Fig. 4.11). In principle, 
this feature may have significant practical applications for all the processes mentioned 
in the Introduction. 



CHAPTER 5 

DEPOSIT GROWTH FOR FINITE- VOLUME PARTICLES 
IN CIRCULAR EVAPORATING DROPS 

In the preceding parts of this work we considered how mass of the contact-line deposit 
grows with time and how it depends on such geometrical characteristics of the drop as 
its radius (circular drops) or its opening angle and distance from the vertex (pointed 
drops). However, we never attempted to describe the geometrical characteristics of 
the contact-line deposit itself, e.g. the width and the height of the deposit ring in the 
case of circular geometry. At the same time, there is solid experimental data (21) 
on various geometrical characteristics of the ring and their dependence on time, the 
initial concentration of the solute, and the geometry of the drop. In this chapter we 
will provide a simple model that addresses this lack of theoretical understanding of the 
geometrical properties of the contact-line deposit and provides analytical results that 
compare favorably to the experimental data. This model is as universal (in its range 
of validity) as all the consideration above and provides the description that depends 
only on the geometry of the problem. This model is based on the assumption that 
particles have finite volume and they simply cannot pack denser than some volume 
fraction, e.g. fraction of close packing. It turns out that this model is sufficient to 
explain most of the observed phenomena. 

The notion that the profile of the deposit could be found by the simple assumption 
that the solute becomes immobilized when the volume fraction reaches a threshold 
was originally suggested by Todd Dupont (39). First efforts to create a model were 
conducted by Robert Deegan (21; 22) who formulated most of the physical assump- 
tions and wrote them down mathematically. Here we present the entire problem, 
including its full formulation and its analytical and numerical solutions (not reported 
in the literature previously). 

5.1 Model, assumptions, and geometry 

In this chapter we will restrict our attention to the circular geometry of the drop only, 
as this geometry is the easiest to deal with mathematically and the most important 
practically. We will continue to use cylindrical coordinates (r, 0, z) with the origin 
in the center of the drop. This will allow us to employ some results obtained for the 
circular drops earlier in this work. 

As we already saw, it is necessary to assume that the contact angle is small in 
order to be able to obtain any analytical results in a closed form. This case is also 
most important practically, as virtually always <C 1 in experimental realizations, 
including experiments of Ref. (21). Thus, we do not lose any generality assuming 
that contact angle is small. This will be the only small parameter in this chapter. It 
is clear that a drop with a small contact angle is necessarily thin, i.e. its maximal 
height is much smaller than its radius. Therefore, it makes good sense to continue 
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considering vertically averaged quantities, like the vertically averaged velocity v of 
Eq. (2.2) that we used in the three preceding chapters. 

Our model pictures the drop as a two-component system (the components are 
"fluid" and "solute"), which has two "phases": the "liquid phase" in the middle of 
the drop and the "deposit phase" near the contact line. Both components are present 
in both phases, and the difference between the phases lies in the concentration of the 
solute. In the deposit phase, volume fraction of solute p is high (e.g. comparable to 
close packing fraction or 1) and fixed in both space and time. Thus, p is just a constant 
number, one can think of it as close packing fraction or 1. (The case of p — 1 may 
seem to be special as there is no water in the deposit phase; however, as will be seen 
below, for small initial concentrations of the solute this case leads to exactly the same 
main order results.) In the liquid phase, volume fraction of solute x varies in space 
and changes with time, and it is relatively small compared to p. The initial volume 
fraction Xi — x(0) is constant throughout the drop; at later moments the solute gets 
redistributed due to the hydrodynamic flows, and the concentration becomes different 
in different parts of the liquid phase. Volume fraction of fluid is then (1 — p) in the 
deposit phase and (1 — x) m the liquid phase. Note that we do not require Xi ^ P m 
this section, although we do assume Xi < V- We should also emphasize that we do not 
presume there is any real "phase difference" between the so-called phases: one phase 
is just defined as having the maximal reachable solute fraction p (solute cannot move 
in this phase) while the other phase is characterized by lower solute fraction x ( m 
this phase solute can move and hence concentration can change in time and space). 
Besides this difference, phases are essentially identical. The idea that the solute loses 
its mobility when its concentration exceeds some threshold was suggested by Todd 
Dupont (39). 

Since the drop is thin and since we employ vertically averaged velocity (i.e. the 
problem is two-dimensional in certain sense), it is natural to assume the boundary 
between the phases is vertical. Thus, the particles get stacked uniformly at all heights 
when being brought to the phase boundary by the hydrodynamic flow (which does 
not depend on vertical coordinate z). This boundary can be pictured as a vertical 
wall at some radius R(t) from the center of the drop, and this wall propagates from 
the contact line (located at Ri = R(0)) to the center of the drop. Fig. 5.1 illustrates 
the mutual location of the phases, and Fig. 5.2 schematically shows the time evolution 
of the drying process and growth of the deposit phase. 

The evaporation rate depends only on the overall shape of the drop, and evap- 
oration occurs in the same fashion from both phases. We assume that evaporation 
is not influenced by any motion of the solute inside the drop, and the necessary 
amount of fluid can always be supplied to the regions of highest evaporation near 
the contact line. Physically, high evaporation near the edge is what brings the solute 
to the contact line, and we assume that presence of the deposit does not obstruct 
the motion of fluid (Fig. 5.3). Since the drop is thin and the contact angle is small, 
we can use expression (3.20) for the evaporation rate, which we used for the case of 
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Figure 5.1: Mutual location of two "phases" in the drying drop: L is the liquid phase, 
and D is the deposit phase. 




Figure 5.2: Time evolution of the deposit phase growth: side view (left) and top 
view (right). Only the deposit phase is shown. Thickness of the ring is exaggerated 
compared to the typical experimental results. 
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Figure 5.3: Presence of particles in the deposit does not obstruct fluid evaporation at 
the edge of the drop. All the necessary fluid is supplied, and it is this motion of the 
fluid that brings the particles to the deposit phase. Also shown schematically is the 
fact that the boundary between the phases is vertical and the particles get stacked 
at full height between the substrate and the free surface of the drop. 



zero-volume particles, but with a slight modification — now, the radius of the drop 
is Ri (while R(t) stands for the radius of the inner boundary of the deposit phase, so 
that R l = R(0)): 

J(r) = l D ^-noo) (51) 



Here D the diffusion constant, and n s and noo are the saturated and ambient vapor 
densities, respectively. The real situation may be different from the assumed above 
when p is large or comparable to 1, and the edge of the area where evaporation occurs 
may be located near the boundary of the phases instead of the contact line. However, 
for small initial concentrations of the solute, the main order result will be insensitive 
to the exact location of the singularity of the evaporation rate: whether it is located 
at the contact line or near the boundary of the phases. We will further comment on 
this case of "dry deposit" when we obtain the full system of equations. 

The geometry of the problem is shown in Fig. 5.4. The radius of the drop is Ri, 
the radius of the phase boundary is R(t), and R(0) = Ri. The height of the phase 
boundary is H(t), and the initial condition is H(0) = 0. We conveniently split the 
height of the free surface of the liquid phase into the sum of H(t) and h(r,t). Since 
H is independent of r, function h(r) satisfies the Young-Laplace equation (2.11). As 
is easy to see, the solution to this equation remains the same as in Chapter 2, and 
hence the shape of the upper part (above the dashed line in Fig. 5.4) is just a spherical 
cap. Thus, we can use expression (3.1) obtained for h(r,t) earlier, which reduces to 
Eq. (3.3) in limit of small contact angles: 

f{ 2 (t) — r 2 



Here 6{t) is the angle between the liquid-air interface and the substrate at phase 
boundary (an equivalent of the contact angle). Note that we do not assume that 6{t) 



53 



0i 



h(r,t) 




Ri 



R(t) 



Figure 5.4: Geometry of the problem. Vertical scale is exaggerated in order to see 
the details, typically H <C and h <ti R{. 



and h(r, t) are necessarily positive at all times: both can be negative at later drying 
stages, and the shape of the liquid-air interface may be concave. Both convex and 
concave solutions for h(r,t) are consistent with Eq. (2.11); the right hand side of this 
equation can have either sign. By definition, both 9{t) and h(r,t) are positive when 
the surface is convex (and therefore they are positive at the beginning of the drying 
process) and negative when the surface is concave. The initial value of 9{t) coincides 
with the initial contact angle 9{ = 0(0). 

Clearly, there are three unknown functions of time in this geometry: 9(t), R(t) 
and H(t). These three time dependences are to be determined in the subsequent 
sections of this chapter. However, these quantities are not independent of each other. 
If we assume that solute particles fill up the entire space between the substrate and 
the liquid-air interface when being brought to the phase boundary, it is easy to see 
that the three geometrical functions are related by the constraint 

"4- = ~<> d 4- (5-3) 
dt dt v 1 

Physically, this assumption means that the angles between the liquid-air interface 
and the substrate are identical on both sides of the phase boundary (9 = \dH/dR\), 
and hence h{r) and its first derivative are continuous past the phase boundary. This 
condition was first introduced by Robert Deegan (21; 22). Thus, there are actually 
only two independent functions of time, 9(t) and R(t). In the case of zero- volume 
particles we had only one independent function of time — contact angle 9(t). 

The geometrical definitions above allow us to find volume of each of the two 
phases. Volume of the liquid phase is simply 

fR(t) fR 3 9 R 2 H\ 

V L = 2vr jf (h(r, t) + H(t)) rdr = 2tt ^— + ^-J . (5.4) 

Taking into account relation (5.3), an infinitesimal variation of this volume can be 
related to the infinitesimal variations of 9 and R: 




dV L = ——d - + 2nHRdR. (5.5) 
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The first term is responsible for the motion of the liquid-air interface, and the second 
term corresponds to the shift of the phase boundary. It is also straightforward to 
obtain an expression for the differential of the volume of the deposit phase, which has 
only the term related to the inward motion of the phase boundary: 

dV D = —2nHRdR. (5.6) 

We will use the last two expressions in the following section. We will also adopt the 
notation that subscripts L and D refer to the liquid and deposit phases, respectively. 

Having formulated assumptions and having stated the model and the geometry, 
we are now in position to write the governing equations. 



5.2 Main equations: conservation of mass 

If one were asked to express the essence of the entire theory in one sentence, this 
sentence would be: "It is all about conservation of mass." Indeed, as we will see 
by the end of this section, all three governing equations obtained here represent 
conservation of mass (or volume) in one form or another. 

We start from the global conservation of mass in the drop. Since the amount of 
solute within the drop does not change, the change of the entire drop volume is equal 
to the change of the amount of fluid only. This fluid gets evaporated from the surface 
of the drop, and the total amount evaporated per unit of time from the entire surface 
of the drop equals to the total change of fluid volume in the drop. Hence the total 
change of drop volume is equal to volume of fluid evaporated from the surface: 

dV ^' = dvF Lf < 5J > 

By convention, superscripts F and S refer to the fluid and the solute components, 
respectively (while subscripts L and D continue to denote phases). Now, the total 
change of drop volume is the sum of volume changes of each phase: 

dV\ tot = dV L + dV D = ^d ^ , (5.8) 

where dV^ and dVjj were found in the preceding section [Eqs. (5.5) and (5.6)]. On 
the other hand, the volume of fluid evaporated from the surface of the drop during 
time dt can be determined from 



<1] 1 = - P J ^J^^2.rdr « -2n P ^ rdr, (5.9) 

surf J P V J P 



dt 



where, as in the rest of this work, p is the density of the fluid and we neglected the 
gradient of h(r) with respect to unity (which is always legitimate for thin drops). 
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Simple integration with J(r) of Eq. (5.1) yields 

4P(n fl -noo)i^ 
surf n 



dV F 



Combining Eqs. (5.7), (5.8), and (5.10), we finally obtain the first main differential 
equation of this chapter: 

R4 d ( 9 \ _ 16D(n s - nopjRj 

dt\R) up 1 ' ' 

This equation represents the global conservation of mass in the drop and relates time 
dependencies of 6{t) and R(t). 

Our next equation represents the local conservation of mass and is analogous to 
Eq. (2.3). Now, one has to take into account that there are two components in the 
liquid phase and write a separate equation for each component. Since solute particles 
are carried along by the flow [a free particle of an appropriate size reaches the speed of 
the flow in about 50 ns in water under normal conditions (40)], the velocities of each 
component should be identical at each point within the liquid phase. Conservation 
of fluid can be written in exact analogy with Eq. (2.3): 

V • [(1 - X )(h + H)v] + l + dt[(l- X )(h + H)] = 0, (5.12) 

P 

where we took into account that the surface height is now (h + H) and neglected 
its gradient with respect to unity in the second term as we did everywhere in this 
work. Similar equation can be written for the conservation of solute, but without the 
evaporation term: 

V-[x(h + H)v} + d t [ X (h + H)} = 0. (5.13) 

In the last two equations x is the volume fraction of the solute at a given point within 
the liquid phase. Adding the two equations and employing the linearity of differential 
operations, we obtain a direct analog of Eq. (2.3): 

V- [(h + H)v] + - + dt(h + H) = 0. (5.14) 

P 

This equation could have been obtained if we considered only one component with vol- 
ume fraction 1 in the liquid phase, and this equivalence should not be of any surprise: 
we explicitly assumed that the solute particles move in exactly the same fashion as 
the fluid does, and hence any differentiation between the two (from the point of view 
of the conservation of volume) is completely lost. Note that if evaporation were too 
intensive, this equivalence would not hold, as there might be an insufficient amount 
of fluid coming into a volume element and the solution could get completely dry (i.e. 
only the solute component would be left). We implicitly assume this is not the case 
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for our liquid phase where the solute fraction is relatively small and the evaporation 
is not too strong. 

In circular geometry Eq. (5.14) can now be resolved with respect to the flow 
velocity, as was done in Eq. (3.22): 

^ *) = - 7U ^ W) /; g + S t h + S tH ) rir. (5.15) 

Straightforward integration with J(r) of Eq. (5.1), h(r,t) of Eq. (5.2), and dH/dt of 
Eq. (5.3) and employment of Eq. (5.11) for d{6 / R)/dt yield 



v r (r,t) = 



2D(n s - rioo) Ri 



7rp 



ii0 



(5.16) 



This is an analog to Eq. (3.23); however, in this expression the velocity does not 
diverge at r — R, but rather reaches a finite value (which is quite natural as the 
evaporation rate does not diverge at r = R). This equation is a direct consequence 
of the local conservation of mass. 

As we did in the preceding chapters, we label r^(t) the initial location of the 
solute particles that reach the phase boundary (and become part of the deposit ring) 
at time t. This function is monotonically decreasing, as the solute particles from the 
outer areas of the drop reach the deposit phase sooner than the particles from the 
inner areas do. The derivative of this function is related to the velocity found in the 
preceding paragraph in a very simple fashion: 



dr; . . 

- = -v r ( ri ,t). 



(5.17) 



We explained the origin of this relation in Chapter 2. Thus, this equation and the 
result for the flow velocity of the preceding paragraph yield the second principal 
equation of this section: 



drj 
dt 



2D(n s 



R0 



1 



(a)' 



(5.18) 



+ H 



This equation relates ri(t) to time dependencies of the geometrical parameters of the 
drop [0(t), R(t), and H(t)]. 

The volume of the solute in the deposit phase v£ at time t is equal to the volume 
of the solute located outside the circle of radius r^t) at time (since all the solute 
between r^(t) and Ri becomes part of the deposit by time t). The latter volume can 
be found similarly to Eqs. (2.15) and (3.25) (which are written for the mass instead 
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of the volume) and hence 




/ h{r, 0) 2-nrdr = V s 1 



( 



Ri 



) 



(5.19) 



where V = nxiRfOi/4: is the total volume of the solute in the drop. On the other 
hand, the volume of the solute in the deposit phase is just the constant fraction p of 
the volume of the entire deposit phase: 



Equating the right-hand sides of these two equations, taking the time derivatives of 
both sides, and making use of the already determined dVp of Eq. (5.6), we obtain 
the third principal equation of this section: 



As the other two principal equations (5.11) and (5.18) are, this equation is a conse- 
quence of the conservation of mass. While Eq. (5.11) comes from the global conser- 
vation of the fluid, this equation represents the global conservation of the solute. 

Thus, we have four unknown functions of time: 9(t), R(t), H(t), and r^(t), and 
four independent differential equations for these functions: Eqs. (5.3), (5.11), (5.18), 
and (5.21). In reality we need only three of these functions: 9(t), R(t), and H(t); 
however, there is no simple way to eliminate r^t) from the full system of equations 
and reduce the number of equations. Having solved this system of the four differential 
equations, we will be able to fully characterize the dimensions of the deposit phase 
and describe the evolution of the deposit ring. The following sections are devoted to 
the details of this solution. 

Here we will only comment on how this system of equations would change if we 
were to consider the case of the completely dry solute p — 1. In this case there 
would be no evaporation from the surface of the deposit phase, and the effective edge 
of the evaporating area would be somewhere in the vicinity of the phase boundary. 
Assuming the same one-over-the-square-root divergence of the evaporation rate at 
r = R instead of r = R^ [which mathematically means substitution of R in place of 
Ri in Eq. (5.1)] and conducting a derivation along the lines of this section, one can 
obtain a very similar system of four differential equations. These equations would be 
different from Eqs. (5.3), (5.11), (5.18), and (5.21) in only two minor details. First, 
Eqs. (5.11) and (5.18) would lose all indices % at all occasions of R4 (i.e. one should 
substitute R for all Ri in both equations). Second, p should be set to 1 in Eq. (5.21). 
Apart from these details, the two systems would be identical. As we will see in the 
following section, this difference between the two systems is not important in the 



V§=pV D . 



(5.20) 




(5.21) 
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main order in a small parameter introduced below, and thus this "dry-solute" case 
does not require any special treatment contrary to the intuitive prudence. 



5.3 Analytical results in the limit of small initial 
concentrations of the solute 

So far we have not introduced any small parameters in this problem other than the 
initial contact angle 0j <C 1. In particular, equations (5.3), (5.11), (5.18), and (5.21) 
were obtained without assuming any relation between p and Xi other than the non- 
restrictive condition Xi < P- I n order to find the analytical solution to this system, we 
will have to assume that Xi <^ P- Then, we will solve the same system of differential 
equations numerically for an arbitrary relation between Xi an d p. 

Assumption Xi ^ P physically means that the solute concentration in the liquid 
phase is small — it is much smaller than the concentration of close packing or any 
other comparable number of the order of 1. This is the case for most practical realiza- 
tions of the ring deposits in experiments and observations: the solute concentration 
rarely exceeds 10% of volume, and in most cases it is far lower. If the volume fraction 
of the solute in the drop is small, then the solute volume is also small compared to 
the volume of the entire drop. Hence, the deposit phase, which consists mostly of 
the solute and whose volume is of the order of magnitude of the volume of the entire 
solute in the drop, must have small volume compared to the volume of the entire 
drop. Thus, if the initial volume fraction of the solute in the drop is small, then 
the dimensions of the deposit ring must be small compared to the corresponding di- 
mensions of the entire drop (i.e. the ring width must be much smaller than the drop 
radius and the ring height must be much smaller than the drop height). 

Let us now introduce parameter e that is small when Xi/P is small. However, we 
do not fix its functional dependence on Xi/P f° r the moment: 

< 1, (5.22) 

where / is an arbitrary increasing function of its argument. Then we postulate that 
the ring width is proportional to this parameter: 

R(t) =Ri[l- eW(t)] , (5.23) 

where W(t) is an arbitrary dimensionless function and we explicitly introduced the 
dimensionality via R{. Obviously, W(0) = 0. So far we have not done anything 
but writing mathematically that the ring width is small whenever the initial volume 
fraction of the solute is small (i.e. the statement of the preceding paragraph). Next, 
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we introduce a dimensionless variable for the angle 9(t): 

0(t) = 9i6(t), (5.24) 

where both 6(t) and Q{ are small everywhere in this chapter, while the newly intro- 
duced function 6{t) is arbitrary (in particular, 9(0) = 1). From the definitions above 
and from the geometrical constraint (5.3), it is straightforward to conclude that for 
small initial concentrations of the solute the height of the ring H(t) must be linear in 
small parameters e and Q{ and directly proportional to the only dimensional scale Rf 

H{t) = e6iRiH(t), (5.25) 

where H(t) is a dimensionless function of time (H(0) = 0). Equation (5.3) fixes the 
relation between this function and the previously introduced functions W(t) and 9(t): 

f=^. (5.26) 
dt dt y J 

Finally, we introduce the last dimensionless variable in place of the fourth unknown 
function rj(t): 

V(t) = l-{^\ (5.27) 

with the initial condition V(0) = 0. So far, we simply introduced four new dimension- 
less variables in place of the four original unknown functions of time and explicitly 
separated dependence of these unknown functions on small parameters e and 0j. As 
a final step of this procedure, we introduce the dimensionless time r: 

r = f , (5.28) 
where t j is a combination of system parameters with the dimensionality of time: 

*f = 16D(n s - noo y (5 ' 29) 

In the limit Xi/P ~^ this combination represents the time at which all the solute 
reached the deposit phase; for finite Xi/p it does not have so simple interpretation. 

Substitution of all the definitions for the dimensionless variables of the preceding 
paragraph into the original system of equations (5.3), (5.11), (5.18), and (5.21) and 
retention of only the leading and the first correctional terms in e lead to the following 
simplified system of equations: 
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dti ^~dW 
dr dr 



~ d6 
3eW— = 
dr 



■1, 



dV 
~dr~ 



Vv - v 2 [1 - AeW (y- 1 - l) 



29V 



1-eW (2V^' 1 -l)] + 4eH 



p dr 



Ae'H^- (l - eW) . 
dr v ' 



(5.31) 



(5.32) 



(5.33) 



These differential equations are still coupled, but the coupling is simpler than in 
the original system. As is apparent from the last equation, parameter e 2 must be 
proportional to Xi/P- Since the separation of the ring width into e and W in Eq. (5.23) 
is absolutely arbitrary, parameter e is defined up to a constant multiplicative factor. 
Therefore, we set this factor in such a way that e 2 is equal to Xi/P : 



(5.34) 



This fixes the function / from the original definition (5.22). 

Next, we expand all four unknown functions of r in small parameter e: 

6 = 9 + e0i + • • • , (5.35) 

V = V + eVi + ---, (5.36) 

H = H + eHi + ■ ■ ■ , (5.37) 

W = W + eWi + ■ ■ ■ , (5.38) 

and keep only the main-order terms after substitution into system (5.30)-(5.33). The 
resulting main-order system of differential equations becomes very simple: 



dH ~ dW 
= #0" 



dr 



dr 



dVo 
dr 



d9o 

dr " ' 

Jv~ -v 2 



20 o Vb 



dr 



dr 



(5.39) 



(5.40) 



(5.41) 



(5.42) 



Clearly, equations decouple: the second equation can be solved with respect to 6o(r) 
independently of all the others, then the third equation can be solved with respect to 
Vq(t) independently of the first and the fourth, and finally the first and the fourth 
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equations can be solved together as well (since dW^/dr can be eliminated from the 
two and the result can be solved with respect to Hq(t)). Thus, it is a matter of 
technical effort and time to obtain the following solution to the system of equations 
above with the appropriate initial conditions: 



o (r) = 1 - r, 



1-r 



,3/4' 



2/3 



Hq(t) 



Wq{t) 



1 



B 



(1-r) 3 / 4 



7 4 

3' 3 



(1 -/)3/4]V3 



(5.43) 
(5.44) 
(5.45) 

(5.46) 



/0 8# (t') (1-r') 1 / 4 

Here B(a, 6) = Jox 0-1 ^ ~~ x)^ 1 dx is the complete beta-function, B z (a,6) = 
Jo i a_1 (l — is the incomplete beta-function (a > 0, b > 0, and < z < 1), 

and the integral in the last equation cannot be expressed in terms of the standard 
elementary or special functions. In a similar fashion, systems of equations of the first 
and higher orders in e can be written, however, they cannot be resolved analytically 
as easily as the main-order system above. 

A system of equations similar to our system (5.30)-(5.33) was presented by Robert 
Deegan in works (21; 22). Small concentrations are assumed in those works, but no 
general system of equations similar to our Eqs. (5.3), (5.11), (5.18), and (5.21) is 
written before introducing the small parameter and expanding in it. The equations 
in that system mix terms of the main and the correctional order in concentration 
arbitrarily. Upon examination, that system does include all the main-order terms of 
our system; however, the first-order corrections are incomplete, and most corrections 
are missing. Moreover, no analysis quantifying how small the neglected and the re- 
tained "small corrections" are is made in Refs. (21; 22), and therefore the importance 
of different terms is difficult to infer from that system. Equations are not written 
as a single system of four differential equations for four unknown functions; instead, 
equations for two variables are solved before even writing the other two equations, 
thereby undermining the equality of all the four functions of time that must be de- 
termined simultaneously and making it impossible to compare the neglected terms. 
Finally, no analytical solution to the system of equations is provided in Refs. (21; 22); 
the author tackles the system numerically. The numerical results are not presented 
explicitly either; they are only used to convert the data between the experimental 
graphs. All these deficiencies are rectified in our approach above. The analytical 
solution (5.43)-(5.46) is reported here for the first time. Higher-order terms can also 
be constructed (with sufficient labor invested), and the procedure can be conducted 
up to an arbitrary order. 
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Time (sec) 

Figure 5.5: Mass of a drying drop as a function of time. Experimental results, after 
Refs. (21; 22). The line running through the data is a linear fit. (Courtesy Robert 
Deegan.) 



How do our results (5.43)-(5.46) translate into the original variables? The first 
two of them [Eqs. (5.43) and (5.44)] reproduce the earlier results obtained for the zero- 
volume deposit. In terms of the original (dimensional) variables Eq. (5.43) represents 
the linear decrease of the angle between the liquid-air interface and the substrate with 
time: 

9(t) =0^1-^. (5.47) 

Taking into account the definition of t-e [Eq. (5.29)], one can see that this is a direct 
analog of Eq. (3.16) for the contact angle in the zero-volume-particle case. So, angle 
9 in the finite- volume case depends on time in exactly the same fashion as the contact 
angle in the zero-volume case does. This expression also provides an interpretation 
of tf\ it is the time at which the free surface of the liquid phase becomes flat. Before 
time tj this surface is convex, after tj it becomes concave and bows inward (until 
it touches the substrate). Thus, tj- is not the total drying time, but rather the time 
at which 9 becomes zero. In the limit Xi/P or equivalently e — ► the height 
of the deposit is going to zero and the two times are the same. For finite values of 
both parameters the total drying time is longer than the time at which the liquid-air 
interface becomes flat. Eq. (5.47) has been verified in experiments (21; 22) where the 
mass of the drop was measured as a function of time (Fig. 5.5). Since the mass of 
the drop is directly proportional to 9 these results confirm the linearity of 9(t) during 
most of the drying process. 

The second equation (5.44) also has a direct analog in the zero-volume case. In 
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Figure 5.6: Fraction of the solute in the deposit phase V^/V as a function of time 
[Eq. (5.49)]. 



the original variables it can be rewritten as 



3/4 



'nity 



3/2 



(5.48) 



which is identical to Eq. (3.24). Clearly, = when t = tf. Thus, tj can also be 
interpreted as the time at which all the solute particles become part of the deposit. 
The same interpretation can be obtained differently: according to Eqs. (5.19), (5.27), 



and (5.44), the fraction of the solute particles in the deposit phase V^/V b is simply 



Yd 
v s 



VA 




4/3 



(5.49) 



(plotted in Fig. 5.6). This fraction becomes 1 when t — tf, and thus all the solute 
particles reach the deposit by time tf. So far, the results of this finite- volume model 
coincide with the results of the zero-volume case considered earlier. 

However, the third and the fourth equations [Eqs. (5.45)-(5.46)] provide com- 
pletely new results. In the dimensional variables they yield the height of the phase 
boundary H and the width of the deposit ring W = Ri — R, respectively: 



H{t) = MOiRiHo ( f 
V P Vf, 



(5.50) 
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Figure 5.7: Height of the phase boundary H(t) in units of OjRi^Xi/p as a function 
of time [Eq. (5.50)]. 



W(t) 



'-RiW 




(5.51) 



where functions Hq(t) and Wq(t) are given by Eqs. (5.45) and (5.46) and plotted in 
Figs. 5.7 and 5.8. (Note that these expressions represent only the leading term in the 
expansions of the full functions for height and width in small parameter e = \jxilv-) 
These results provide the sought dependence of the geometrical characteristics of the 
deposit ring on all the physical parameters of interest: on the initial geometry of the 
drying drop [R{ and 9{), on the initial concentration of the solute (xi), and on the time 
elapsed since the beginning of the drying process (t). If the time is considered as a 
parameter, they also provide all the necessary information to obtain the geometrical 
profile of the deposit (i.e. the dependence of height on width), which we plot in 
Fig. 5.9. Note that the vertical scale of this plot is highly expanded compared to the 
horizontal scale since there is an extra factor of 9{ <C 1 in the expression for height; 
in the actual scale the height is multiplied by Q{ (in addition to the same scaling 
parameters as in the width) and hence is much smaller than it appears in Fig. 5.9. 

The time dependence of H and W plotted in Figs. 5.7 and 5.8 deserves a brief 
discussion. It is straightforward to obtain the asymptotics of Hq(t) and Wq(t) for 
small values of the argument and for values around 1 (i.e. for the early and the late 
drying stages). At early times, both the height and the width scale with the drying 
time as a power law with exponent 2/3: 



(3r) 2 /3 




[l + 0( 




(5.52) 
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Figure 5.8: Width of the deposit ring W(t) in units of RiyXi/p as a function of time 
[Eq. (5.51)]. 




Width W/(eR.) 

Figure 5.9: Deposit ring profile: height vs. width. The vertical scale is different from 
the horizontal scale by a factor of Q{ <C I. 
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(5.53) 



Thus, at early times H « ^W^, which can also be deduced directly from Eq. (5.3) 
without obtaining the complete solution above. At the end of the drying process, 
the height and the width approach finite values (which, apart from the dimensional 
scales, are universal, i.e. simply numbers), and do so as power laws of (tf — t) with 
two different exponents: 
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where Hq(1) and Wo(l) are just numbers: 
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(5.54) 
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Clearly, dH/dW — 0^(1 — r) and hence vanishes when r — > 1. This fact can also be 
observed in flattening of the graph in Fig. 5.9 at late times. 

Dependence of the height and the width on the radius of the drop R^, while 
intuitively obvious (since Ri is the only scale in this problem with the dimensionality 
of the length), has been verified in experiments (21; 22). A linear fit has been obtained 
for the dependence of the ring width on the radius (Fig. 5.10), which agrees with our 
findings. 

Comparison to the experimental data for the dependence on the initial concentra- 
tion of the solute is slightly less trivial. Our results predict that both height H and 

1/2 

width W scale with the initial concentration as (at least, in the leading order 
for small concentrations). However, experimental results by Deegan (21; 22) show 
a different exponent of Xi- Exponent of Xi was determined to be 0.78 ± 0.10 and 
0.86 ± 0.10 for two different particles sizes (Fig. 5.11). Why is the difference? The 
answer lies in the fact that the width measured in experiments (21; 22) is not the 
full width of the ring at the end of the drying process, but rather the width of the 
ring at depinning. Depinning is a process of detachment of the liquid phase from the 
deposit ring (Fig. 5.12). This detachment was observed experimentally in colloidal 
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Figure 5.10: Ring width vs. drop radius. Experimental results, after Refs. (21; 22). 
The line running through the data is a linear fit. (Courtesy Robert Deegan.) 



suspensions but has not been explained in full theoretically yet. An important ob- 
servation, however, is that the depinning time (i.e. the time at which the detachment 
occurs and the ring stops growing) depends on the initial concentration of the solute. 
This dependence was also measured by Deegan (Fig. 5.13). The resulting exponent 
was determined to be 0.26 ±0.08. Thus, the width of the ring at depinning scales 
with the initial concentration of the solute Xi as 

W d oc X }/ 2 W (fy oc £*Wo (x?- 26± °- 08 ) . (5-58) 

where is the depinning time (t^ftf oc ^-9- 26=l= 0- 08 ). As is apparent from Fig. 5.13, 
the typical values of the depinning time are of the order of 0.4-0.8 tj. In this time 
range, function Wq is virtually linear (Fig. 5.8). Therefore, the dependence of on 
Xi has the overall exponent of the order of 0.5 + (0.26 ± 0.08) = 0.76 ± 0.08. It is 
now clear that both experimental results 0.78 ± 0.10 and 0.86 ± 0.10 fall within the 
range of the experimental uncertainty of the approximate predicted value 0.76 ±0.08. 
All in all, it turns out the theoretical dependence of the ring width on the initial 
concentration agrees with the experimental results very well. 

Note that Robert Deegan (21; 22) did not report direct measurements of the 

1. While the full explanation is yet to be developed, the naive reason for depinning seems rela- 
tively straightforward. The pinning force depends only on the materials involved and is relatively 
insensitive to the value of the contact angle. At the same time, the depinning force is simply the 
surface tension, which is directed along the liquid-air interface and which increases as the contact 
angle decreases (since only the horizontal component of this force is important). Thus, the relatively 
constant pinning force cannot compensate for the increasing depinning force of the surface tension, 
and after the contact angle decreases past some threshold, the depinning force wins and causes 
detachment of the liquid phase from the deposit. 
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Figure 5.11: Ring width normalized by the drop radius vs. initial concentration of the 
solute for two different particle sizes. Experimental results, after Refs. (21; 22). The 
two data sets are offset by a factor of 5 to avoid mixing of the data points related to 
the different particle sizes. The lines running through the data are linear fits in the 
double-logarithmic scale, which upon conversion to the linear scale yield power laws 
with exponents 0.78 ± 0.10 and 0.86 ± 0.10. (Courtesy Robert Deegan.) 




Figure 5.12: A photographic sequence demonstrating a depinning event. Experimen- 
tal results, after Refs. (21; 22). The view is from above, and the solid white band in 
the lower part of the frame is the ring; the rest of the drop is above the ring. The 
time between the first and the last frames is approximately 6 s; the major axis of the 
hole is approximately 150 /zm. (Courtesy Robert Deegan.) 
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Figure 5.13: Depinning time normalized by the extrapolated drying time vs. initial 
concentration of the solute. Experimental results, after Refs. (21; 22). The line 
running through the data is a linear fit in the double-logarithmic scale, which upon 
conversion to the linear scale yields a power law with exponent 0.26 ±0.08. (Courtesy 
Robert Deegan.) 



height of the deposit. He calculated the height from the data in hand, and thus direct 
comparison to the experimental data is not available for the height [and comparison 
between the plots of Refs. (21; 22) and our predictions depends on the details of 
Deegan's calculation and its inherent assumptions]. 

The square-root dependence of the height and the width on the concentration is in 
good agreement with general physical expectations. Indeed, the volume of the deposit 
ring is roughly proportional to the product of the height and the width. On the other 
hand, the height is of the same order of magnitude as the width since the ratio of 
the two is of the order of &i (which is a constant). Thus, both the height and the 
width scale approximately as a square root of the ring volume. Finally, the volume 
of the deposit ring is proportional to the initial volume fraction of the solute: the 
more solute is present initially, the larger the volume of the deposit ring is at the end 
(since the ring is comprised of only the initial amount of the solute). Therefore, both 
the height and the width must scale as a square root of the initial volume fraction. 
It is rewarding that the results of our complex calculation match the results of this 
simple physical argument. 

Thus, the complete analytical solution to our model is available in the limit Xi/p ~~ * 
0, and this solution compares quite well to the experimental results. Since the main- 
order solution in Xi/p is perfectly adequate, the difference between the original system 
of equations and the one for the "completely dry" case is not important: in the main 
order in Xi/P the results are identical for both cases (since one case in different from 
the other only by presence of R instead of in a few places in the main equations, 
and this difference is of the correctional order in Xi/p)- 
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5.4 Numerical results for arbitrary initial concentrations of 

the solute 

Apart from approaching the original system of equations (5.3), (5.11), (5.18), and 
(5.21) analytically, we also solve the same system numerically. During this numerical 
procedure we do not presume that Xi/p is small, nor do we expand any quantities or 
equations in e or any other small parameters. Our main purpose is to reproduce the 
results of the preceding section and to determine the range of validity of our analytical 
asymptotics. 

Typical values of Xi/P m m ost experimental realizations are of the order of 0.001- 
0.01, and thus, only concentrations below approximately 0.1 are of practical interest. 
Thus, we will mostly concentrate on this range of practical importance while de- 
scribing the results despite the fact the numerical procedure can be (and have been) 
conducted for any ratio Xi/P- ^ i s to be noted that in the case of Xi comparable to p 
our model is not expected to produce any sensible results, as the entire separation of 
the drop into the two phases (the liquid phase and the deposit phase) is based on the 
assumption that the mobility of the solute is qualitatively different in the two regions. 
When Xi is comparable to p the two phases are physically indistinguishable, while the 
model still assumes they are different. Therefore, it should be of absolutely no sur- 
prise if any results for Xi ~ P are unphysical or unrealistic. When the assumptions of 
a theory are explicitly violated, its results are not to be taken seriously in that range 
of parameter values. For completeness, we do provide results for x% — Pi however, in 
no way we attempt to claim these numerical results describe the real physical state 
of the system at this value of the parameter. As we mentioned above, only the values 
of Xi/P around 0.1 and below are presented as our final numerical results. The case 
Xi = P is provided only as an illustration of the general trend, without any attempt to 
draw any conclusions from this unphysical (within our model) value of the parameter. 

We present our numerical results for the same quantities (and in the same order) 
as in our analytical results (5.47), (5.49), (5.50), and (5.51). Since for arbitrary Xi/P 
time tf is not exactly the total drying time, there is a question of where (at what 
time) to terminate the numerical curves. By convention, we terminate all the curves 
(except Xi/p — 1) m all the graphs at value of t/tj when all the solute reaches the 
deposit phase. In our model, it turns out that this time approximately coincides with 
the time the center-point of the liquid-air interface touches the substrate. For all the 
initial concentrations (except Xi = P)i the time the center-point touches the substrate 
was numerically found to be within 0.1% of the time the very last solute particles reach 
the deposit phase. Thus, within our model, the moment the center-point touches the 
substrate and the moment the last solute particles reach the deposit ring are about 
the same, and the curves are terminated at exactly this time. Of course, in reality a 
small fraction of the solute stays in the liquid phase as long as the liquid phase exists, 
and so the moment the last solute particles reach the deposit phase must be after the 
moment the center-point touches the substrate; however, the amount of the solute 



71 




Figure 5.14: Numerical results: Dependence of angle 9 between the liquid-air interface 
and the substrate at the phase boundary on time t. Different curves correspond to 
different initial concentrations of the solute; values of parameter Xi/p are shown at 
each curve. The analytical result in limit Xi/p ~ > is also provided. 



remaining in the liquid phase at touchdown is insignificant, and thus practically all 
the deposit has already formed by that time (99.9% of all the solute is already in the 
ring). 

Numerical results for angle 9 as a function of time are shown in Fig. 5.14. All 
curves (except Xi/P = 1) behave almost linearly (as expected), however, the slope 
increases with the concentration: formation of the ring in the drops with more so- 
lute finishes faster [in the relative scale of tj of Eq. (5.29)]. The end of each curve 
demonstrates the value of the angle 9t at the moment the liquid-air interface touches 
the substrate. (The analytical expression for this angle is 9t = —2H/R for a thin 
drop.) The absolute value of this angle increases with concentration, which is quite 
natural since the height of the ring grows as a square root of the concentration while 
the radius of the liquid phase does not change substantially for small concentrations. 
Clearly, the numerical results converge to the analytical curve when Xi/p ~^ 0- 

Growth of the volume fraction of the solute in the deposit phase with time is shown 
in Fig. 5.15 for various solute concentrations. This graph reconfirms the observation of 
the preceding paragraph that the solute transfer happens faster (in units of the time 
scale tf) for denser colloidal suspensions. All curves are terminated when volume 
fraction Vp/V S becomes equal to 1. The apparent termination of the curve for 
Xi/p = 0.1 earlier than that is an artifact of the plotting software, which discarded 
the very last point of this data set when creating a plot. Presumably, this plot and 
the preceding one should hold true independently of the geometrical details of the 
solute accumulation in the ring (which cannot be expected from the following plots 
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Figure 5.15: Numerical results: Dependence of the volume fraction of the solute in 
the deposit phase Vp/V S on time t. Different curves correspond to different initial 
concentrations of the solute; values of parameter Xi/P are shown at each curve. The 
analytical result in limit Xi/P ~ ¥ ls a l so provided. 



for the ring height and width). 

The next two graphs represent the numerical results for the geometrical character- 
istics of the deposit ring as functions of time: the height is shown in Fig. 5.16, while 
the width is in Fig. 5.17. The ring profile, i.e. the dependence of height on width (ob- 
tained by elimination of time from the two results), is also provided in Fig. 5.18. As 
the graphs depict, the ring becomes wider and lower at higher initial concentrations 
of the solute. Since the volume of the ring is roughly proportional to the product 
of the height and the width, the decrease in height must be of the same percentage 
magnitude as the increase in width. This can be qualitatively observed in the graphs. 

As a final piece of the numerical results, we create a double-logarithmic plot of 
the dependence of the height and the width on the initial concentration of the solute 
(Fig. 5.19). The predicted square-root dependence on the initial concentration is seen 
to hold true for volume fractions up to approximately lCP 1 / 2 ^ for the height and up to 
approximately lCP^/ 2 p for the width. Deviations for higher volume fractions are due 
to the increasing role of the correctional terms in e compared to the main-order terms 
in expansions (5.37) and (5.38). In this graph, as in all the results of this section, it 
is clear that our main-order analytical results provide an adequate description of all 
the functional dependencies in the range of the initial concentrations of experimental 
importance (0.001-0.01). 

In general, the numerical results of this section complement and reinforce the 
analytical results of the preceding section, and the good agreement between the results 
of the two independent solutions provides a crosscheck of both methods. 
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Figure 5.16: Numerical results: Dependence of the height of the phase boundary H 
on time t. Different curves correspond to different initial concentrations of the solute; 
values of parameter Xi/P are shown at each curve. The analytical result in limit 
Xi/p — > is also provided. 
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0.0001 
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Figure 5.17: Numerical results: Dependence of the width of the deposit ring W on 
time t. Different curves correspond to different initial concentrations of the solute; 
values of parameter Xi/p are shown at each curve. The analytical result in limit 
Xi/P ~ > is also provided. 
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Figure 5.18: Numerical results: Ring profile. Dependence of the height of the phase 
boundary H on the width of the deposit ring W . Different curves correspond to 
different initial concentrations of the solute; values of parameter Xi/P are shown at 
each curve. The analytical result in limit Xi/P ~ > * s & l so provided. 




Figure 5.19: Numerical results: Log-log plot of the dependence of the height of the 
phase boundary H and the width of the deposit ring W on the initial volume fraction 
of the solute. "Lg" stands for the logarithm of base 10, and thus the numbers along 
the horizontal axis represent the order of magnitude of the ratio Xi/P- The main-order 
analytical results H oc yXi/p an d W oc yXi/p are a l so provided for comparison. 
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5.5 Discussion 

The specific graphs of Figs. 5.14-5.18, both analytical and numerical, represent the 
original contribution of this chapter and may be determined experimentally giving 
validation to the model and our theory. Measurements of the profiles in Fig. 5.18 
should be particularly easy to conduct (as there is no time dependence involved) and 
may confirm or refute the predicted robustness and universality of the deposition 
profiles. 

One may notice that the curves in Fig. 5.18 end at some positive (non-zero) 
height. This indicates the solute is exhausted before the profile curves had a chance 
to return to the substrate, and the final shape of the deposit ring must have a vertical 
wall at its inner end. We believe this is an artifact of our model that is inherently 
two-dimensional when flows inside the drop are concerned. Thus, we used the depth- 
averaged velocity (2.2) throughout this work, and we also explicitly assumed the 
phase boundary is vertical and the particles get stacked uniformly at all heights (i.e. 
the vertical distribution of the solute was assumed homogeneous). This is equivalent 
to assuming that vertical mixing is complete. This assumption is quite important, 
and the results are expected to get modified if the true velocity profile (2.7) (or any 
other three-dimensional distribution) is used instead of the depth-averaged velocity. 
We expect that if a three-dimensional model were built and the dependence on z 
were taken into account for all quantities then the discontinuous wall of the phase 
boundary would get smoothened and the height would continuously return to zero. 
A question remains whether such a model would be solvable analytically. 

Our model relies on the assumption that the mobility of the solute is different in 
the so-called liquid and deposit phases. Essentially, we assume that the mobility is 
in the deposit phase and 1 in the liquid phase. This assumption, while artificial in 
its nature, seems relatively reasonable when applied to this system. Indeed, in the 
physical situations near the close packing, the loss of mobility typically occurs over 
a quite narrow range of the concentration values, and hence our assumption should 
work satisfactorily when the difference between \i an d P is in the orders of magnitude. 
The higher the initial concentration is and the closer the two values are, the worse this 
assumption holds true and the more artificial the difference between the two phases 
is. We observe this in our numerical data: the numerical results for Xi = P are quite 
unphysical, while the low-concentration results seem to be consistent and coherent 
with expectations. Thus, the validity of any model based on this separation of the 
mobility scales decreases for higher initial concentrations of the solute. 

The model assumes that the free-surface slope between the liquid and the deposit 
phases is continuous. In fact, assumption (5.3) expressing this continuity is one of 
the four basic equations of this chapter. This assumption seems quite natural as well. 
Indeed, if the liquid is present on both sides of the phase boundary, the change in the 
slope of its free surface would cost extra energy from the extra curvature at the phase 
boundary, since the liquid-air interface possesses effective elasticity. Presence of this 
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extra energy (or the extra pressure) at the location of the phase boundary is not 
justified by any physical reasons under the conditions of this chapter where processes 
are slow and the surface is in equilibrium. In equilibrium the surface shape must 
have constant curvature past the phase boundary since the entire separation into the 
two phases is quite artificial as discussed above. Presence of the particles below the 
liquid-air interface does not influence the surface tension, and thus the liquid surface 
(and its slope) should be continuous at the phase boundary. On the other hand, if 
the density of the particles matches the density of the liquid (which was the case 
in the experiments), nothing prevents the particles from filling up the entire space 
between the substrate and the liquid-air interface, thus providing the growth of the 
upper edge of the deposit phase along the liquid-air interface. This is particularly 
true for the thin drops discussed in this chapter (^ <C 1 was assumed throughout) 
where the vertical mixing is intensive, where the free surface is nearly horizontal, 
and where the problem is essentially two-dimensional. However, the equality of the 
slopes on both sides of the phase boundary does not seem inevitable, and one may 
think of the situations when it does get violated. One example might be the late 
drying times when the deposit growth is very fast (see Fig. 5.15) and hence the 
deposition may occur in some non-regular manner inconsistent with the slow-process 
description of this chapter. Other examples may be related to the gravity (slightly 
unequal densities of the particles and the liquid) or the convection. In any case, this 
assumption can possibly be checked experimentally, and if condition (5.3) is found 
violated, an equivalent constraint dependent on the details of the deposit-growth 
mechanism must be constructed in place of Eq. (5.3) in order to relate functions 6{t), 
R(t), and H(t). 

Another inherent assumption of our model is related to the evaporation rate J(r). 
As we discussed at the beginning of this chapter, we assumed that presence of the 
solute inside the drop does not affect the evaporation from its surface. This is gener- 
ally true when the evaporation is not too fast and the deposit phase is not too thick 
and not too concentrated. When these conditions are not obeyed, presence of a thick 
or concentrated layer of the solute on the way of the liquid moving from the phase 
boundary to the contact-line divergence of the evaporation rate may create a strong 
viscous force. This viscous force would prevent the necessary amount of the fluid 
from being supplied to the intensive-evaporation region near the contact line. Gen- 
erally, we assumed throughout this work that the viscous stresses are not important, 
and showed in Chapter 2 and the Appendix that this assumption is valid whenever 
v <C o~/3r). In the deposit phase, the velocity is large due to the proximity to the 
contact-line divergence of the evaporation rate, and the effective viscosity is large due 
to the high concentration of the solute. Thus, this condition may get violated and the 
viscosity may become important in the deposit phase, slowing down the supply of the 
liquid and ultimately making the deposit dry. Obviously, this affects the evaporation 
rate, and the functional form of the evaporation profile changes. Simple assumption 
that the evaporation rate stays of the same functional form, but with the divergence 
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at the phase boundary (at R) instead of the contact line (at Rj), was shown above 
not to affect our main-order results. Thus, our results appear to be relatively insen- 
sitive to the exact location of this divergence within the (narrow) deposit phase. (In 
reality the evaporation edge would be somewhere between the contact line and the 
phase boundary, i.e. the real situation is intermediate between the two considered.) 
However, the deposit could modify the evaporation rate J(r) in other ways. When 
there is a dry deposit ring just outside the liquid phase, the entire functional form of 
J may change, and the Laplace equation for an equivalent electrostatic problem must 
be solved anew with the additional boundary conditions responsible for the presence 
of the dry solute rim and the modified evaporation at the edge. As we already saw in 
all the preceding chapters, this is the most complicated part of the problem, and the 
mathematics can become prohibitively complex. Thus, finding the exact form of J 
may be a formidable task. One way around is in creating such evaporating conditions 
that the functional profile is simpler, for instance, J is just a constant. This would 
be more difficult to control experimentally, but would be much easier to treat analyt- 
ically. The unavailability of the exact analytical form for J seems to be the biggest 
open question in this class of problems (the same difficulty was also encountered for 
the pointed-drop solution of Chapter 4). 

The equilibrium surface shape of the liquid phase is a spherical cap (5.2). This is 
a rigorous result proved in the preceding chapters and valid during most of the drying 
process. However, when h(0,t) becomes negative and exceeds H(t) in absolute value 
(i.e. when the center-point touches the substrate), the surface shape is no longer 
spherical. Moreover, a new element of the contact line is introduced in the center 
of the drop in addition to the original contact line at the perimeter, and the entire 
evaporation profile gets modified in addition to the modified surface shape, thus 
influencing all the other quantities. Our treatment does not account for the small 
fraction of the drying process occurring after this touchdown (which is a change 
in topology of the surface shape, and thus requires a separate treatment after it 
happened). First of all, the amount of liquid remaining in the drop at this moment 
is of the order of e compared to the original volume and hence would not modify 
our main-order analytical results. Second, as our numerical calculations show, at 
touchdown practically all the solute is already in the deposit phase, and the remaining 
amount of the solute in the liquid phase is insignificant. Thus, within our model, the 
remainder of the drying process cannot modify the deposit ring substantially, and 
hence this neglect of the late-time regime seems well justified. Experimentally, the 
inner part of the deposit ring is different from our prediction (which is a vertical wall) 
and appears to have a spread shelf. Presence of this tail in the deposit distribution 
can be caused by several features absent in our model. Its inherent two-dimensionality 
may be one of these shortcomings (as discussed above); the account for the dynamical 
processes occurring after the deposit phase has already been formed (e.g. avalanches 
of the inner wall) may be another missing feature. Absence of the account for the late- 
time regime after the center-point touches the substrate may be among these reasons 
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influencing the final distribution of the deposit as well. A more detailed account for 
the effects of this late-time regime might be required in the future. 



CHAPTER 6 
CONCLUSIONS 



The major results presented in this work are related to two problems. One of them 
is the problem of the solute transfer and the deposit growth in angular evaporating 
drops, and the other is the problem of determining the geometrical characteristics of 
the deposition patterns in circular evaporating drops. 

For the angular drops, we provide the full solution to the problem of the deposit 
accumulation at the contact line. We determine all the necessary physical quantities 
(the surface shape, the evaporation rate, the flow field, and the mass of the deposit) 
and fully describe the deposition process near the vertex of the angle. Not surprisingly, 
we find that the angular drops are rich in singularities that govern all these physical 
quantities and affect the deposition. All the quantities scale as power laws of the 
distance to the vertex and the elapsed drying time with exponents dependent only 
on the opening angle of the drop. These exponents are universal and do not depend 
on any free or fitting parameters, they depend only on the drop geometry. This 
dependence on the geometry is useful for creating controlled deposition patterns and 
is a step towards complete understanding of the evaporative deposition for arbitrary 
contact-line geometries. 

For the problem on the geometrical characteristics of the deposition patterns, we 
provide a model that accounts for the finite dimensions of the deposit on the basis of 
the assumption that the solute particles occupy finite volume and hence these dimen- 
sions are of the steric origin. Within this model, we find the analytical solution for 
small initial concentrations of the solute and the numerical solution for arbitrary ini- 
tial concentrations of the solute. We demonstrate the agreement between our results 
and the experimental data and show that the observed dependence of the deposit 
dimensions on the experimental parameters (the initial concentration of the solute, 
the initial geometry of the drop, and the time elapsed from the beginning of the 
drying process) can be attributed mainly to the finite dimensions of the solute parti- 
cles. These results are also universal and important for understanding the deposition 
process and controlling the pattern formation. 

The entire work provides all our results and derivations pertaining to the evap- 
orative deposition and obtained to date in a single volume, which may serve as a 
reference guide to anyone interested in the theory of the evaporative deposition. We 
hope this reference volume will simplify future theoretical and experimental efforts 
towards the complete understanding of these ubiquitous phenomena. 
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APPENDIX A 
DOMINANCE OF THE CAPILLARY FORCES OVER 

VISCOUS STRESS, AND LEGITIMACY OF 
EMPLOYMENT OF THE EQUILIBRIUM SURFACE 

SHAPE 

The purpose of this section is to demonstrate that for sufficiently slow flows one can 
employ the equilibrium surface shape for finding the pressure and the velocity fields 
instead of having to solve for all the dynamical variables simultaneously. We will also 
quantify how slow "sufficiently slow flows" are. 

We start from the equation of the mechanical equilibrium of the interface (2.10), 
where we approximate the doubled mean curvature 2K with V h. This approximation 
holds true because the liquid-air interface is nearly horizontal for thin drops and small 
contact angles (this can be proved rigorously for both geometries of interest, see the 
corresponding sections in the main text), and the other terms of the functional K[h] 
are unimportant. Substitution of 

p= -<TV 2 h + Patm (A.l) 

into the Darcy's law (2.8) yields 

v = v*h 2 V{V 2 h), (A.2) 

where v* = a/3r). Upon further substitution into the conservation of mass (2.3), we 
obtain 

V • \v*h 3 V(V 2 h)} +- + d t h = 0, (A.3) 
J P 

which, together with Eq. (A.2), constitutes the full system of equations for finding 
h(r, 0, t) and v(r, 0, t). 

Now, for water under normal conditions, rj — 1 mPa-s and a = 72 mN/m. Hence, 
the velocity scale v* is of the order of 

v* = ^- pa 24 m/s. (A.4) 

Obviously, this is a huge value compared to the characteristic velocities encountered 
in usual drying process. Therefore, one can develop a systematic series expansion in 
small parameter e = v/v* (where v is some characteristic value of velocity, say, 1 or 
10 //m/s): 

h = ho + ehi + --- + e n h n + ■■■, (A.5) 
v = v + evi + • • • + e n w n + ■■■, (A.6) 
and keep only the Hq and vq terms at the end in order to describe the process up to 
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the main order in e = v/v* . A similar expansion can also be constructed for pressure: 

p = po + m + ■ ■ ■ + £ n Pn + ■■■, (A.7) 

where pi are related to hi by equation (A.l): 

PO = -o-V 2 h +p atm , pi = -aV 2 hi, etc. (A.8) 

Physically, condition v <C v * is equivalent to the statement that the viscous stress is 
negligible and that the capillary forces dominate. Let us understand what Hq and vq 
physically correspond to. 

Plugging the expansions for h and v into the system (A.2)-(A.3), one obtains a 
set of terms for each power of e, starting from and up. Equating terms of the 
main order in e yields the following two equations 

fyjV(V 2 /i ) = and V • [fyjV(V 2 /i )] = 0, (A.9) 

which both can be satisfied if and only if V 2 /io is a function of time only. Writing it 

as 

w2l PO ~ Patm 1 /a in\ 

vhQ = — — = "w (A - 10) 

we immediately identify this equation with the statement of spatial constancy of the 
mean curvature of the interface, which describes the equilibrium surface shape at any 
given moment of time t (i.e. we obtained equation (2.11) with the desired properties 
of pq). Thus, ho is indeed the equilibrium surface shape, and so is h (up to the 
corrections of the order of v/v*). 

Repeating the same procedure for the terms of the next order in e, we arrive at 
another two equations: 

v =^§V(V 2 /ii), (A.ll) 



/i§V(V 2 /ii)] +- p +d t h = 0, (A.12) 

which can be seen to be equivalent to the set of equations (2.12)-(2.13) upon identi- 
fication tjj = vV 2 h\ = —epi/3i]. Knowing the equilibrium surface shape ho, one can 
solve the second equation above with respect to the reduced pressure ip, and then 
obtain velocity vq by differentiating the result according to the first equation. Thus, 
up to the corrections of the order of v/v*, one can first find the equilibrium surface 
shape h(r, <fi) at any given moment of time, and then determine the pressure and the 
flow fields for this fixed functional form of h, as was asserted in the main text. 
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